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PREFACE 

This  is  the  forty-seventh  in  a  scries  of  I’eports  gi  owing  out  of  the  study  ot 
radar  cross  sections  at  The  Radiation  Laboratory  of  The  University  of  Michigan. 
Titles  of  the  reports  already  published  or  presently  in  process  of  publication  are 
listed  on  the  preceding  pages. 

When  the  study  was  first  begun,  the  primary  aim  was  to  show  that  rada ' 
cross  sections  can  be  determined  theoretically,  the  results  being  in  good  agreement 
with  e^eriment.  It  is  believed  that  by  and  large  this  aim  lias  been  achieved. 

In  continuing  this  study,  the  objective  is  lo  cetermine  means  for  computing 
the  radar  cross  section  of  objects  in  a  variety  of  doierent  environments.  This  has 
led  to  an  extension  oi  the  investigation  to  include  not  only  the  standard  boundary- 
value  prddems.  but  also  such  topics  as  die  emission  and  propagation  of  electro¬ 
magnetic  and  acous*^ic  waves,  and  phenomena  connected  with  Ionized  media. 

Associated  witti  Jie  theoretical  work  is  an  eiqierimental  program  which 
embrac-*^  (a)  measurement  of  antennas  and  radar  scatterers  in  order  to  verify  data 
determined  theoretically;  (b)  iuvestigetioa  of  antenna  behavior  and  cross  section 
problems  not  amenable  to  theoretical  solution;  (c)  problems  associated  with  the 
design  and  development  of  microwave  absorbers;  and  (d)  low  and  hi^  density 
ionization  phenomena. 

R.E.  Hiatt 
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I 

INTRODUCTION 

This  is  the  first  of  a  series  of  reports  aimed  at  summarizing  tne  available 
information  about  the  sc::ltering  properties  of  selected  bodies  of  simple  shape. 
Perhaps  the  sittipl>->ir  of  all  shapes  is  tlie  sphere,  and  it  is  probable  that  more  has 
been  written  about  this  one  body  than  about  all  other  bodies  put  ^ether.  To  detail 
ail  of  the  results  in  one  repo.r^  is  therefore  impossible,  and  iu  seekhig  to  summarize 
them  so  as  to  provide  an  intelligible  account,  an  author  is  compelled  to  restrict  him¬ 
self  to  those  theories  and  those  methods  of  sciutiun  which  he  feels  are  most  signi¬ 
ficant. 

In  taking  as  the  subject  of  this  first  report  the  diffraction  of  electromagnetic 
energy  by  the  sphere,  our  object  is  to  gathev  together  in  one  place  some  of  the  more 

j 

useful  forms  of  solution,  both  exact  and  approximate,  giving  also  a  brief  account  of 
the  methods  of  derivation.  Wherever  possible  references  are  given  to  tabulations  of 
tiie  functions  and  series  involved,  particularly  ocro'.cstion  with  the  standard  Mie 
solution.  I 

I 

Section  II  ir-  devoted  tc  the  Mie  solution  and  since  this  is  the  starting  point 
for  most  ot  vhe  other  approaches,  a  detailed  descripliou  is  given.  Certain  special 
applications  are  discussed,  and  references  are  given  to  computationr  based  on  thr.  ! 


Mie  series. 
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For  sufficiently  low  frequencies  an  alternative  representation  of  the  solution 
is  possible  in  which  the  field  components  are  e^qianded  in  ascending  positive  integral 
po'vers  of  ka,  where  k  is  the  wave  number  and  a  is  the  radius  of  the  sphere.  The 
corresponding  expansion  for  the  far  field  amplitude  is  the  s  v  laiied  Rayleigh  series, 
and  this  is  described  in  Section  Ill.  Two  derivations  are  given:  in  the  first  of 
these  the  series  is  obtained  by  expanding  the  various  terms  in  the  Afie  solution,  but 
in  the  second  the  low  frequency  expansion  is  obtained  directl}*  without  any  eiqiiicit 
reference  to  the  Mie  result. 

Section  IV  is  concerned  with  the  high  frequency  scattering  behavior  and  the 

approach  which  is  adopted  is  based  on  the  Watson  t-«msform.  In  recent  years  the 

Watson  transform  technique  has  teen  generalized  to  an  extent  which  permits  the 

asymptotic  solution  of  a  large  class  of  diffraction  problems,  and  the  general  method 

stems  from  the  fact  that  locaiiv  all  convex  iiodies  with  radii  of  curvature  much 

j 

larger  than  the  wavelenf  'h  are  similar  to  a  sphere  of  radius  equal  to  that  of  the 
conve;"  '-•ody  in  the  direction  of  enei^  flow.  This  local  analysis  led  Fock  |ip4^  to 
construct  certain  universal  functions  which  have  been  computed  and  tabulated  by 
Logan  ^195^  .  Since  this  material  dues  not  appear  in  any  standard  refei  ence,  a 
relatively  detailed  eiqxisition  is  given. 

In  the  final  section  the  physical  optics  approach  is  considered  insofar  es  it 


applies  to  tne  sphere  problem.  The  approximate  expressions  for  the  current 


THE  UNIVERSITY  OF  MICHIGAN 

3G48-1-T 


disti'ibution  and  for  the  Tar  field  are  compared  with  the  exact  expressions  dt  rived 
from  the  Mie  series,  and  a  numerical  comparison  for  ka  ~  1C  is  presented. 

Aithou^  the  previous  methods  have  covered  the  region  of  'applicaliility*  of  physical 
optics,  it  was  felt  that  this  approximate  but  well-known  technique  should  be  included 
both  for  historical  eucrest  and  because  of  the  ease  with  which  rough  and  read)- 
answers  can  be  obtained  thereby.  For  a  general  and  more  critical  ejqmsition  of 
the  physical  optics  method,  the  reader  is  referred  to  Baker  and  Copson  [l95o]. 
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THE  EXACT  SOLUTION 


This  seciiun  is  devoted  to  the  exact  6''iution  of  the  problem  of  scattering  of  a 
(dans  electromagnetic  wave  by  a  sphere.  A  brief  account  of  the  derivation  is  in¬ 
cluded  since  almost  all  subsequent  computations  and  approximations  rdy  to  some 
extent  on  this  exact  result.  In  addition  to  the  homogeneous  sphere  and  the  important 
limiting  case  of  perfect  conductivity,  results  for  two  concentric  siAeres  are  also 
presented.  The  sim^difications  stemming  from  the  "far  field"  assumptluii  arc  also 
discussed.  A  guide  to  computed  results  available  in  the  literature  is  presented  and 
some  representative  curves  are  included. 

2. 1  The  Mie  Series  For  the  Sphere 


Tlie  first  exact  solution  for  the  scattering  of  a  plane  wave  by  a  homogeneous 
sphere  is  usually  attributed  to  Mie  ^90^  although  much  work  was  done  before  then. 
Thompson  flSSsJ  treated  tlie  perfectly  conducting  sphere  with  equal  rigor,  and  in 
his  exhaustive  work  on  the  sphere  Logan  |i95^  gives  precedence  to  Clebsch  pSCsJ  . 
Nevertheless  the  series  solution  for  the  st^ere  in  terms  of  spherical  wave  functions 
is  usually  referred  to  as  the  Mie  series  and  this  general  usage  w'lll  be  employed 
here.  Descriptions  of  the  solution  abound  in  the  literature,  the  most  popular,  per- 
fhaps,  being  that  given  by  Stratton  |l94lj  and  it  is  his  presentation  on  which  me 

I 

■present  account  is  based. 
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The  problem  is  that  of  determining  the  electric  and  magnetic  field  vectors, 
E=E*  +  E®,  andH=H*+H® 

(where  i  and  s  denote  incident  and  scattered  respectively),  external  to  a  homo¬ 
geneous  sphere  of  radius  a,  permcabilit}'#4i,  permittivity  tj,  and  ccxiductivity  Sj 
in  the  presence  of  an  incident  or  primary  field  given  by 

e'-eSc-*", 

"  °  (2-1)* 

„i  „  /V  -ikz 
H  = y  e 

A  rectangular  Cartesian  coordinate  system  (x,  y.  7)  has  been  employed  in  which 

eqns  (2-1)  describe  aplane  wave  travellingin  the  direction  ofthc  negative  z-axis  with  its 

1  k 

electric  vector  confined  to  the  x  direction.  H  =  Y  E  .  where  Y=  -z  - - is  the 

00  Z  uUq 

intrinsic  admittance  cf  free  space;  k  is  the  profagation  constant  of  liic  medium  in 
which  the  sphere  is  imbedded,  wliich  medium  is  assumed  homogeneous,  isotropic, 
and  a  perfect  rliclectrit  and  is  here  taken  as  free  space.  In  terms  of  the  permittivi^ 


and  p  .•ni“ability, 


1  ! - -  2sr 


where  A  is  the  wavelength.  M.k.s.unit3  are  employed  and  the  harmonic  time  factor 
c  has  been  suppressed.  The  restriction  to  free  space  a  trivial  one  because  in 
a  mediu;i*.  characterized  by  e  and  p  different  from  their  free  space  values,  a  pro* 
Mgation  constant  k  may  be  defined  as 


An  underlined  symbol  denotes  a  vector  and  a  caret  denotes  a  unit  vector. 
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k  =  uycjr  .  (2-3) 

Similarly,  if  the  conductivity  s  is  non-zero,  the  propagation  constant  can  be  taken  as 

k  =  4i(e  +  ^)  .  (2-4) 

It  is  convenient  to  have  the  center  ol  the  sphere  coincide  with  the  origin  of  the  coor¬ 
dinate  system.  This  detracts  none  of  the  generality  and  permits  the  use  of  splicrical 
polar  coordinates  (r  6,  0)  where 


x  =  rsinOcos^,  y  =  rsin6sin0,  andz  =  rcos0, 
in  terms  of  which  the  surface  of  the  sphere  is  sitr.pl;  i  ”  a  (see  Figure  2-1  ). 
The  free  space,  source  free,  Ma.xwell  equations  are 

I 

i 

S  bB 

V.n-c,  — =0, 


V-  H  =  V*  E=  0  . 


After  sufHiressing  the  harmonic  time  I'arialion  these  equatimis  require  that  all  field 

quantities  e.xtcrior  to  the  sphere  he  solutions  of  the  vector  wave  equation. 

V„(V-F)  =  k2F  .  (2-7) 

i  s  i  s 

••  .«re  F  can  be  E  ,  E  ,  H  or  H  .  Interior  to  the  si^erc,  E  and  H  must  satisfy 

VA(V*F)-kiF.  (2-U) 

where  k,  is  the  propagation  constant  for  the  material  comprising  the  Sfdiciv. 
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Tbe  boundary  conditions  are  really  continuity  conditions  at  the  surface  oftfae  sphere,  i.e., 

=  rAEl 


taJe'+E®) 


r=a 


r=a 


(2-9) 


rA(H*+H®)| 


H  +H  )  =rAH| 

lr=a'‘'  lr=a“ 

General  solutions  ihe  vector  tvavc  equation  can  be  generated  by  vector 
operations  on  the  solutions  of  the  scalar  aavc  equation 

{V*  +  k*)t>  =  0,  (2-10) 

in  the  following  way.  If  c  is  a  solutitw  of  fcqn(j£-10)  th«i  the  three  vectors 


L 

il  =Va(Ev) 


(2-11) 


are  orthogonal  solutions  of  eqn  (2-7).  These  are  known  as  Hansen's  vector  wave 
hincUoiu..  having  been  proposed  by  Hansen  [1935.1936,1937]  in  his  woric  on  radiation 
from  antennas.  They  are  discussed  more  fully  by  S'jarton  ^19^l)and  Senior  jlSB^  . 

Since  field  quanuties  are  required  by  Maxwell's  equati<ms  to  be  solenoidal, 

[or  divergence  free,  the  fact  that 

V-L=^d  =  -k?v>¥=0  .  (2-12) 


I 


shows  that  only  the  M  and  N  vectors  can  be  involved  in  their  representation. 
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hhc  following  form: 


*  rib  . 


n(frrl)  ,,  .  frt  ctx  ! 

<2-16? 

«ntrc  tnc  pnircs  indicate  dilTerentiation  with  respect  to  kr.  the  subscripts  c  aw3 
.  are  short  for  "enen-  and  ’odd-  respectively  and  refer  to  the  0  dependence  of  the 
chariicierisuc  solution  c  .  and  the  superscript  d)  denotes  the  radial  fanctim;  used. 

The  «iK:.-.ript  <2)  wil!  be  used  to  denote  the  funr-n.  rt^^ainedif  j^(kr>  is  rc|daced 
by  h^<kr)  a.Hl  Ci)  wU  be  used  if  ykr)  is  replaced  by  j^(^r). 

Since  the  field  quantities  are  sdenoidal  the  most  general  e.vpression  for  the 
scattered  electric  field  is 

OP  CO 

y  ^  lA  -.n  v<2)  . 

~  e  4  ^  (2~n. 

m-o  n-o  j^mn  ^mn  oHin^mn 
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v/htTc  the  coefficients  A  and  H  inr:A  .  only  Uie  pr^uagation  constants  k 

®nin  ®inn 

o  o 

and  kj  and  the  sj^crc  radius,  a  .  From  ^taxweli’r  '3  {2-6},  it  is  seen  then 

that 

rvs  tw\ 

'  .  (2-18) 


s 

ia=o 


.  (2) 

f  (B  m'  +A  N' 

^  c  -e  e  “c 

!Pc  mn  mn  mn  mn 
o  c  o  o 


Kliere  the  coefficients  in  this  equation  arc  tLe  same  as  in  eqa{2-ll). 

S:mi!ar!3'.  tne  most  goieral  expressions  for  the  fields  within  the  sphere 
(r<a).  are 

E  =  E  ^  (C  +D  )  .  (2-19) 

and 


m=o  n=o  mn  mn  mn  mn 
o  c  o  o 


^(3),e  n<3))  . 

~  UU'  L _ t  i—t  c  “  e  e  “  e 

*  m=o  n^  mn  mn  mn  ron 
o  o  o  o 


{2-2C) 


l^ere  the  constants  C  and  B  usaiu  ia\x>lve  wily  the  {mopagation  constants 

mn  mn 

*'0 

land  radius. 

die  incidimt  fidd  given  by  eqn  (2-1).  expansions  in  tcm;s  of  vector  wave 
functimis  arc  given  in  Bratton  [f d4lj  as 


_ia 


e‘=  E  )  (-i)"  (m'V  +i  ) 

“  o  n(n+l)  oin  eln 


(2-21) 
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£  H)”  — -  (  N* 


o  n  =  l 


n{n  +  l) 


eln  ola 


(2-22) 


No  terms  corresponding  to  n=0  occ-tr  ^x^causc  (cos  B)  s.0. 

Strai^tfonrard  sultstitutiun  of  cqns  (2-17)  u>  (2-22)  in  the  continuity  rclattons 

(2-3)  no’-v  leads  to  the  following  \'aluos  for  the  unknotrn  coefficients: 

A  =D  -C  =D  -0.  for  all  m  and  n, 
cmn  oii:ii  emn  omn 

A  =B  =C  =D  =0,  for  m  ^  1  and  all  n. 
omn  emn  omn  cmn 

oln  i^n+1)  pj  s^(jija)^k3h^(ka)J  -p^h^{ka)^aj_(^aj 

eln  iKn+i)  ^^|^(kj3)|}^h^(ka^  h^(ka)[kja  y;k,aH 

(2n+l) «, 

-  - -  — — —  '  — '  —  '  ' 

kan(n+l)|^tj^h^(ka}^k.aj^(k|3)j  -pj  j^(lqa)^kah^(kaj  J 

jj  .  _ (-i)”  (2r.-«-l)pi  k|  a _ 

”  a(trt-l)|pj(ka)^h^{ka^k|aj^(kja)j  -#.'^(kia)'J^(l4a)|tah^<ka)J  J 

(2-23) 


where  the  prime  denotes  UiffcrcntisUon  wJUi  respect  to  ka  or  i^a  as  apfTTjpriatc. 

The  situation  is  cofssickrably  simj^cr  in  the  important  case  when  ‘In*  coa- 

•* 

ductivity  of  the  sphere  becomes  infinite  (Im  ki-*0D  ).  The  ctmtinuity  cwiditien 
r^n  (2-9)  is  then  rei^aced  by  the  bmiitdary  condition 
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=  0  . 


sinct-  no  fields  can  exist  aithin  the  s}»«re,  airf  making  use  of  ifce  a^'niptrtic  forms 
of  the  sjxieric&l  Bessel  and  Kankc!  Auctions,  the  cocBicicnta  in  cqns  (2-23)  become 

.  .  j_an  ML. 

*  oin  *  h  (ka) 

n 


hai  (lia)] 

-  n _ J 


aln+l)  |kah^(fcaj 


<2-25] 


C  =  D  =0. 
aln  eln 


2. 2  The  Mie  Series  For  Taro  C^mtceatric  ^ 


The  complicated  proolen  resultisg  iraen  tltc  Sf^re  is  not  hooK^cneous 
i«t  consists  of  a  tonK^raeous  s|**'re  coterwl  «ith  a  honsc^ncoas  layer  at  diffident 
material,  has  bet*  s«ved  by  Aden  airf  Kerter  JiOSlJ  .  The  geometry  is  esscntmily 
the  same  at  pictared  in  Fi}«re  2-1.  e3mept  for  the  mUition  oS  a  ^ifaec  1^-cr  of 
thickness  d,andtliis  is  shown  in  Figure  2-2  wtere  Ur* ,  ?s:ih*ex-axis  and  Incii^t  t 
licld  point  out  !hc  psgo. 

Ouislstcnt  with  tte  notatimi  of  the  prertous  section,  the  iisicr  sphere  of 
tadius  a  will  be  diaract<  rized  by  ki,  €|.  |a*|,  Sj:  the  layer  by  k2,  €3.  p,*  S’- 
ICKi  the  whole  spherical  strucnire  ctf  radius  b  -  a+d  will  be  indtetUed  ir  'ree 
s^ce  ci.aractcrized  by  k,  msd  . 


13  — 


d 


xnc  vi^iiviiivoxix  t. 

3G48-1-T 


FIGURE  2-2: 


In  each  region  the  representations  of  the  field  quantities  are  different.  For 
the  inner  sphere  and  for  free  space,  the  representations  are  similar  to  those  used 


in  section  2. 1,  viz 
r>b 


E=£*+£®,  H=H^+H® 


_i  ,  ,  .,n  2n->-l  ,,il)  ,  ,,(l)v  , ,,  .,n  ?n+l  ,..,(1}  ,Jl). 

E  =£  /  (-i)  (m  +1  N  ,  H  =-iH  >  (-i)  -j--  t<iM  +N„|_), 

-  o^— y  n(n+l)  “oln  “eln  -  ^  n(n+l)  — eln 

n=l  n=i 


£®=E  ;  (A  +B  ).  H®=-iH  >  (B  +A  ) 

-  o  1—^  n-olr.  n~eln  “  o/  n '■eln  n  “oln 


(2-26)+ 


For  conv<'nience  the  coeffici/  nts  are  written  A  ,  B  ,  etc. ,  rather  than  H 

n  n  oln 

B  ,  etc. 
eln 
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E=  E  ^  (C  +D  ),  i!._i  (D  ^,(3)  . 

n-oln  n-ein'*  uwi  Z_  '  n-eln  n-oln^  ' 

while  in  the  layer  r<  b 

E=£  M^f +;^? 

n-cln  n-oln  Si -eln  n  ^In  ’ 

H  =  -i  ^  ^  */2 

^  ^n-eln  "Si-W^n-oln  %-oln  * 


(2-27; 


(2-28) 


I 

where  the  superscripts  on  the  wave  functions  indicate  the  radial  functioic  ’.viiich 

occur.  Tlius  and  defined  in  eqn  (2-16),  '"ntain  the  radial  function 

ykr).  If  this  is  replaced  by  h^(kr),  and  result.  Similarly,  replacing 

i^(kr)  by  jj^(kiv}  yields  and  replacing  j^(kr)  by  jjj(kzr)  yields  and 

and  by  h^(k2r)  gives  M  and  The  continuity  relations  require  r>vE  and 

H  to  be  continuous  a  the  interfaces  r  =  a  and  r  =  b.  This  provides  a  sufficient 

numbe-  of  equations  to  determine  the  unknown  coefficients,  of  'vhich  onjv  A  and 

^  n 

B^,  the  coefficients  of  the  scattered  field,  are  presented  here.  Tlicy  arc 


'l* 

kb  j  (kb)]' 
L  ■’n 

r 

^2 

h  (kb)f  +j 
n  1  1 

kbh  (kb) 

L  n  ^ 

T - 

^2 

B  =  (-i) 
n 


,n+l  2n+l 


jn(kb)fj[kbykb)j' 


(2-29) 


n(n+l)  h  (kb)f7f(kbh  (kb)!'  f 
n  3  n  _  4 
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where 


j  (kja) 

r  -  A--  - 
1 


||^k2bj^^(k2b)j  -|^k2aHk2a^  |{2b  h  ^(k2b)j 'J 


^  ~  [l<2l>l‘^(k2Wj  -h^(k2a)[k2bj^(k2b)]  '|  • 

~"[[Mjn(k2a^  hjl^b)-[k2ah^(k2a,^  jjj(k2b>| 


ij  ikiaM  r  1 

- 1  j  jj(k2b>  h^(k2a)-j^(k2a)  h  Jk2b)j  . 

f  ^  [kiaj^(kiaj  |j^(k2a)[k2bh^(k2bj -h^dtja)  [l^bykjb)]  j 

"  j„(k,a)[[k2byk2bj’  [k2ah^(k2a)j’-[k2aj^(k2a)J[k2bh^(k2bj'|  » 

1 

=  -p-p  [kja  j^^(kia)]  |yk2b)h^0^a)-j^(k2a)h^(k2b)j 


Tr  V  r  V  1 

7iW  |[Mj„(k2a))  h^(k2b)-[l:2ah^(k2aj  yk2b)j  . 


(2-30) 


Scharfman  |^1954j  considered  the  limiting  case  when  the  Inner  sphere  becomes  per¬ 
fectly  conducting  and  the  continuity  condiuon  at  this  interface  is  replaced  by  the 
boundary  condition  eqn  (2-24),  The  expressions  (2-29)  for  the  coefficients  oi  the 
scattered  field  are  still  valid  but  the  f's  defined  by  (2-30)  simplify  as  foilows- 
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1  j  "  ^  jjj(k2a)^k2bh^(k2bj  -h^(k2a)  |^k2b  j^(k2b;j  | 

f„  =  —  {  j  (k2b)h  (k,a)-j  (k2a)h  {kiWh 

2  n  -  n  ^  n  j 

‘l  "  iT  (k'l  |(|^2hjj^(k?b)J[k2ah^(k2aj'-  [k2aj^(k2aj'  [kjok  ^(kzbj'j 

|^hn(k2b)[k2ajn(k2a)]'-yk2b)[k2ah^(k2a)jj  .  (2-31) 

2.  3  The  Mte  Series  For  Two  Concentric  Disjoint  Si’^iercs 

Wben  the*  inner  n'dius  of  the  layer  Is  larger  liian  the  radius  of  the  inner 
sphere,  i.  c.  a  plane  wave  is  incident  upon  a  si^cre  with  two  layers  of  different 

I 

linaterial  upon  it  (sec  Figure  2-3),  the  problem  is  even  more  complicatct!. 

i  r  T 

I  Plonus  jl9GlJ  has  treated  this  proldein,  though  not  in  complete  generality. 

;The  problem  is  sp^eiali/t  J  in  the  following  ways:  1)  the  inner  sphere  of  radius  a 
is  perfev  siv  conducting;  2)  the  two  regions  a  <  r  <  b  and  r  >  c  consist  of  the  same 
material  (here  taken  as  free  space  and  characterized  by  k,  Cq,  and  p^),  and 
3)  the  permeability  of  the  layer  b<  r<  c  is  also  taken  to  be  although  iht  propagatioi 
constant  k2  is  different  from  k.  The  procedure  is  exactly  the  same  as  before.  There 
will  be  three  representations  of  ihe  field  in  the  three  regions,  two  continuity  con- 
Uii.ons  (at  r=b  and  r=c)  and  one  boundary  condition  (at  r=a). 
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FIGURE  2-3 


Thus  for, 
r>n 

E=E*  +  e®,  and  + 

where  the  quantities  arc  exactly  as  defined  in  eqns  (2-2P) 
c>r>h 


E  and  H  are  given  by  eqns  (2-28)  with  ^2=/^ 
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b>r  >  a 


L=E  >  OnM+ZiM+T^N+iN,  , 
o  “  “oln  n  “oln  n”eln  n”eln 

H  = - —  ^  y  +  a  +fl  . 

-  uu  n  “eln  n -eln  n-oin  n  "oln 


<2-32) 


At  thu  interfaces  (r=b  and  r=c),  and  H  must  he  continuous  and 

n  A  •  These  conditions  provide  a  sufficient  number  of  equations  to  deter¬ 

mine  the  unknowTi  coefficients.  Again  only  the  coefficients  of  the  scattered  field. 


A  and  B  are  presented  here: 
n  n 


2n4-l 

j  (kc)f,  + 
n  1 

^kcjn(.ke;j  ig 

n(nH) 

h  (kc)f,  + 
n  1 

kc  h  (kc) 

L  n  J  2 

B  =  -K-i) 
n 


n+1  2nH 


j|/kOf^[kcykc)]’f^ 


(2-33) 


li  (kc)f„  +rkc  h  (kc)l’f . 
n  3  1  n  J  4 


where 


fk2cj^(k2c)j  |h^(ka)^h^(k^)jkbj^(ld>|  -j^(kb)|^k^  h^ik^b^ ^ 

-j^(ka)^h^(k2b)|idj  h^(kbj  -hn(kb)jk2b 
-^k2ch^(k2cj  ^h^(ka)^j^(k2b)  fkbj^(kb]^'-j^(kb)^k2bj^(k2b)|  ^ 

-j^(!a)^j^(k2b}  [kb  h^(kb)J-h^(kb)  [kjb  j^(k  ’ 
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,Uka)f  j^(k2c)f  h^djjb)  [kb h^Ckblj’-h^a*) [k^ 

-h^(kiv)fj^<k2b)fl:bh^(khj’-h^(kb)[k2bi^(J:2bj')| 

-hn(ka)^  jj,(k2c)  [kbj^{kb)j  -j^(kb)[k2bhj[k2b)]’) 

MUk,e)^j^(k2b)[kbj^(kb)]‘-ykb)[kibyk2bj^  . 

h^<>*)[l<2bM!<2b)j-h^(k2b)[kbh^(kb))’) 
-[kacl^Kjcj  ^  ■n^ti<b)[k2bj^(k2b)j -j^(k2b)[kbh^(kb)J^l 
-[k^h„<kaj' [[kac  j^(k,cj  Y(  ^  j  j„Hl^  h^Ckibj^ 

-[k2ch^(k2cj]^(|^Tyi*)[k2bj^(k2bj’-j^(k2b^  . 

V  Mfczb)[kb.ykbj’) 

-fkaj„(ka)j’(h^(kb(k^h^(k,b)JX^)^„^^^ 
-h^<k2c)|^[kah^(kaj|'  (ykb)[>«2bj,<i«2bi|’-('jj'  yk2b)[kbykbj) 
-[kayka)J(h  to(k^i^(k,b)J-(|J .Uk3b(kbh^(id,^)} . 


(2-34) 


When  b=a,  these  expressions  "o  over  to  those  given  in  eqas  (2-31)  fo*-  p  =u 
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2. 4  The  Far  Field  Amplitude 

Of  particular  interest  is  the  far  field  or  far  zone  behavior  of  the  scattered 
field.  Recall  that  the  scattered  field  (exterior  to  the  sphere  and  any  layers)  is 
always  written  as 


i  ^  (A 

o^^  n  — oln  n  -  eln 


(B  ) 

n“eln  n“oln 


(2-35) 


where  the  A  and  B  arc  eiven  bv  eqns  (2-23)  (2-29),  or  (2-33)  depending 

n  n  “  ■ 

on  which  particular  sphere  problem  is  being  considered.  Regardless  of  how  these 
coefficients  are  defined  the  expressiems  for  the  wave  functions  M  and  N  can  be 
simplified  in  the  far  field  of  the  sphere  and  its  layers  (if  any).  Specilicaliy  the 
spherical  HarJeel  functions  contained  in  the  expressions  for  M  and  N  can  be  re¬ 
placed  by  ihe  first  terms  in  their  asymptotic  expansion  for  large  argument  and 
since 

ikr 

h  (kr)~(-i)  ^ -i  —  fkrh  (krl  ,  (2-36) 

n  kr  kr  L  n  J 

the  6  and  0  components  of  both  M  and  N  are  of  eq^al  order  for  giv^i  n.  By  coin- 
I^rison,  tne  radial  component  of  N  is  of  one  hi^er  order.  Consequently,  oruy 
I  tl:e  9  and  0  components  can  appear  in  the  far  field  which  then  has  the  form  of  a 
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It  is  now  con'^onient  to  introduce  the  concept  of  a  scattering  function  rt6, 0). 
Tliis  will  be  defined  by  the  equation 


E®-E  —  . 

-  o  kr 


(2-40) 


valid  for  r»kc*.  whereT  is  a  unit  vector  in  the  direction  of  E‘  and  accordingly 
f(0, 0}  can  be  regarded  as  the  far  Held  amplitude.  From  Maxwell's  equation  we 


then  have 


iKr  .  * 

\-f(O.0)?^r  . 


(2-41) 


and  consequently  the  same  function  describes  both  the  electric  and  magnetic  fidds. 
For  scattering  in  the  backward  direction,  the  scattering  function  will  be  written  as 
f(0),  since  there  is  no  dependence  on  0,  and  cqn  (2-38)  then  gives 


B  )  . 

n 


(2-42) 


Similariy,  for  scattering  in  the  forward  direction. 


f(jr)r.^  i"'^  (A  -iB  )  . 


(2-43) 


(see  eqn  (2-39)  ). 

The  definition  of  f(O,0)  given  in  e^  (2-40)  differs  from  that  usually  ?.do|Med. 
'Ihe  function  f/k  corresponds  to  the  scattering  lunctior.  more  commonly  defined,  but 
ih's  has  the  dls^vantage  of  not  being  dimensionless.  In  clectroma-gnetic  theory 
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(and,  indeed,  in  all  branehes  of  physical  science)  there  sc-errs  to  be  c-vcr\'  advan- 
t^e  ailached  to  using  a  non-ditr.ensional  function,  and  it  is  for  this  reason  that 
the  present  definition  has  been  chosen  in  spite  of  the  fact  that  it  represents  a 
break  from  conwniional  notation.  As  defined  above,  the  funetimi  f  is  indep«ident 
of  r  and  can  be  likened  to  a  polar  diagram  factor.  It  depends  on  angular 
variables  6  and  0  and  on  the  properties  of  the  scattering  body,  and  is  sufficient  to 
specify  the  far  field  in  its  entirety. 

It  is  a  simple  mati>-i‘  to  calculate  the  scattering  cross  section  in  terms  of 
the  function  f.  llie  differential  cross  sc-ct:or.  or  histatic  radar  cross  section  o(6t  0) 
is  defined  by 

Iii|  ' 

and  hence,  from  eqn  (2-4u), 

‘^<^'.0)=  §  •  <2-45) 

An  al'emative  expression  is 

0(0.0)-^  lf(0.0f  .  (2-46) 


«-here  X  is  the  wavelength,  and  the  dimensions  of  a  are  here  made  explicit. 

The  total  scaUering  cross  section  is  related  to  o(0, 0)  by  the  et.uation 

or„=  ^  /o(e,0)d  £2  \2~4D 

where  d  £2  is  an  element  of  solid  an^e,  aix)  by  inserting  eqn(2-45)  we  now  have 


i 
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o^=  jVj(f(e.0^'dn  .  <2-48) 

An  additional  relation  between  and  f  is  provided  by  the  "forward 
scatterinji  theorem".  This  was  first  discovered  in  atomic  theory  and  since  then 
its  electromagnetic  cciui\*alcnt  has  received  a  variety  of  i:>dependent  prottfs  <see  , 
fur  e.xample.  Schiff  ^95^  ,  Jones  |l930  ,  and  de  Ho(^  |l95^  ).  Hie  theorem  is 
merely  an  c.xprcssion  of  conservation  of  energy  and  leads  to  the  equation 

Im.  «»)  .  (2-49) 

where  Im.  denotes  the  imaginaiy  part.  In  atklili.-in  to  the  scattering  function  and 
cross  section  defined  above  there  exist  many  quautiues  in  the  literature  with 
similar  names  but  different  definitions.  This  unfertunate  situation  is  virtually 
uncorrcctalile  at  this  stage  and  the  best  one  can  do  is  exercise  care  in  checking 
}  definitions  and  be  re.^igned  to  the  fact  that  many  existing  results  may  require 

I 

;  rcnormalirati.m  hefore  use.  Some  of  the  more  common  tpiantities  are  presented 
here. 

If  the  scattered  far  field  is  written 

E®=  E„  0+  E  -  ^  .  <2-30) 

~  9  9 

where  E^  and  are  defii  ed  in  eyn  (2-37).  then  uikM:  components  can  be  e.xpressed 
io  the  form 
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ikr 

ikr 

E®  =  -E  sin  0  SJi6)  . 
0  o  kr 


(2-31) 


where  Si(0)  and  %(9)  are  defined  by  referring  to  eqn  (2-37)  and  are  called  the 
complex  amplitudes  of  the  scattered  radiation  for  the  two  poiarizationc. 

The  squares  of  the  absolute  values  of  Sj  and  &  are  called  the  intensities  of 
scattered  radiation  for  the  two  polarizations. 


The  absorption  cross  section  a  and  the  scattering  cross  section  a  are 

2L  ^ 


defined  as 


a 

""a  P. 
1 


P. 


(2-32) 


(2-53) 


where  P  is  the  power  absorbed  by  the  obstacle,  P  the  power  scattered,  and  P, 
a  SI 

the  pawer  incident.  If  no  power  is  absorbed  in  the  obstacle  ami  the  surround!:^ 
medium  is  non-dissipative  (e.  g.  free  space),  then  o^  is  the  same  as  defined 
above. 

The  aim  o  ■  -o  is  known  as  the  extinction  cross  scctism  and  in  cases  where 
a  s 

is  mm  zero  cqn(2-4a),  the  forward  scattering  tbeorem^nust  be  altered  to  read 


o  a  =  ^  Im.  f(s) 
a  s 


(2-34) 
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T A BLE  II-l:  LIST  OF  AVAILABLE  MIE  THEORY  FUNCTIONS 


_  .  Index  of 

Reference  „  ,  ..  H 

Refraction 

Values  of 
ka 

Quantity  Calculated 

Shoulejkin 

1.32 

1,3,00 

Scattering  functions  every  20° 

Blumer  |i92i^  , 

[l93l] 

1.25 

1.33 

1.3 

CO 

0.4,0.S,  1.6, 4.8 

1.0.3 

4 

0.1.0.5.1.3,5.10 

Scattering  functions  every  10° 

Stratton  and 

Houchton  (I93)! 

1.33 

0-40 

Scattering  coefficient 

CasMrsson  [193^ 

|i93.'i 

1.63 

l.riG 

1.50 

0.71-3.16 
(17  values) 

Scatteri^  functions  at  0°,  45°, 
90°,  135  .  and  180° 

Gumpncht,  Sung,  Chiu, 
and  Slieocevich  11952] 

1.33 

6. 8.  i0-35(5) 

Scattering  functions  every  10° 

Gumpricht  and 
Sliepcevich 

1,33 

1.44 

1.20 

20. 30.40.60.80 
100,200,400 
20,80,150 

20.80 

Scattering  coefficient 

KerKer  and  Perlee 

I1953J 

2.00 

1.30-2.80 
(12  values  not  in 
Lowan  tables) 

Scattering  functions  at  90° 

Kerker  and  Cox  ji95^ 

2.00 

3.0-b.C  (11  values  not  Scattering  functions  at  130° 
in  Lowan  tables 

Engelhard  and 

Freiss  TiOStI 

.1.44 

0.4, 1,1.5, 2, 2.5, 3, 

4. 3,8 

Scattering  functions  every  10° 

Paranjpe,  Naik,  and 
Vaidva  09391 

1.33 

4,5,0.7,8,9,10,12, 

20.30 

Scatteringfun*  lions  every  10° 

Ruedy  (l94:ij  {1944] 

1.33 

1/8, 1/4, 3/8, 1/2. 
-3/4.1 

Scottering  coefficient 

(continu{«J  o”  next  page) 
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TABLE  II-l  (continued) 

IndcM  of 

Values  of 

Quantity  Calculated 

Reference 

Refraction 

ka 

Houghton  and 
Chalker  1949 

1.33 

"-24  (33  values) 

Scattering  coefficient 

Lowan  1948 

1,33 

0.5-G.0  (15  values) 

Scattering  functions  every 

1.44 

10°  and  scattering  cotTfficient 

1.50 

2.00 

4.21  -  2.51i 

0.100-1 .00(. 05) 
1.0-3.0(.l) 

Extinction  coefficient 

5.55-  2. 85i 

0.10-1.00(.05) 

1.0-2.0(.l) 

Extinction  coefficient 

8.18  -  1.96i 

0.100-1.00(.025) 

Extinction  coefficient 

3.41  -  1.94i 

0.10-1.0G(.C5) 

Extinction  coefficient  and 

1.0-5,G(.l) 

A  and  B 
n  n 

7,21  -  2,65i 

0.1000-1.000(.025) 

1.00-1.30(.05) 

Extinction  coefficient 

8.90  -  .69i 

0.10-0.30(.01) 

Extinction  coeiicient  and 

0.300-0.430(.005) 

A  and  B 

0,43-0.60(.01) 

n  n 

Riley  1949 

1.486 

0.5-3.0(.)) 

Scattering  functions  every 

10^  and  scattering  coefficients 

Aden  1950 

9.01 -0.  Ji 

0.6-60 

Scattering  functions  at  0° 

00 

0. 6-6.0 

Gumpricni  and 

1.20 

1-6(1) 

and  B 
n  n 

Sliepcevich 

1.40 

8 

1951 

1.50 

10-100(5) 

1.60 

100-200(10) 

200-400(50) 

1.33 

4. 5. 6. 8 

l.-'-4 

10-100(5) 

100-200(10) 

200-400(50) 

(continued  on  next  page) 
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TABLE  n-1  (continued) 


„  ,  Index  of 

Reference  , 

Refraction 

Values  of 
ka 

Qumitity  Calculated 

Kerker, 

Langleben  and 

Gunn  |i93^ 

0.126 

Backscattering by  particle  con- 
sistuig  two  concentric  spheres . 
Inner  sphere  m=l. 75,  outerspherc 
m=8.9-1.5i 

Scharfman  jl95^ 

1.26 

Backscattering  by  lossless  dielec¬ 
tric  coated  perfectly  conducting 
sphere.  Outer  sphere!.  6<m^oo  . 

Hey,  Steivart,  Pinson  oo 
and  Prince  [l95^ 

(K.oDlo 

Backscattering  function  f(0).^ 

The  actual  quantity  tabulated  ia  ^2,  not  f  as  listed  at  head  of  each  column 
(see  Hey  and  Senior  ^95^). 


TABLE  U-2:  PARTIAL  LISTING  OF  TABLES  TO  BE  FOUND  IN  VAN  DE  HIJLST 


Refractive 
Indc.;.  r.' 

25ra/X 

Page 

No. 

Quantitv  Calculated 

CO 

.1(.1)1.6 

1.8-90 

161 

Efficiency  factor 

.8..53,lte 

1.33.1.5.2 

178 

Maxima  and  minima  of  the  extinction 

uui've. 

m  close  to  1 

- 

180 

Extinction  and  absorption  by  partially 
ab.<sorblne  spheres. 

complex 

- 

273-274 

Complex  values  of  m  for  v-hich  computa¬ 
tions  have  been  made. 

3.41-1.941 

7.20-2.651,00 

1.3 

2Ti 

Extinction  coefficient  and  intensity 
functions. 

1.50-in' 

(n*  small) 

.5-7.0 

295 

Extinction  by  spheres. 
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TABLE  n-3:  PARTIAL  LISTING  OF  GRAPHICAL  PRESENTATIONS  OF  DATA  IN 
VAN  PE  HULST  [195“^ 


Refractive 
Index,  m 

2jraA 

Page  No. 

Content 

2 

0-4 

137 

Phase  angle  vs  2valX 

2 

1.55 

1.5 

1.44 

1  33 

1.25 

0-12 

151 

Extinction  curves  of  sphere 

152-153 

Scattering  diagrams 

■BB 

00 

0-3 

162 

Eff  loi-  •  tcy  factors  for  extinction  and  for 
radiation  oressure. 

00 

lA-lO 

163 

Scattering  diagrams. 

1.5 

1.33 

1  +  e 
.93 

0-20 

177 

Extinction  curves  computed  from  Mie's 
formula. 

1.33 

10 

236 

Scattering  diagrams. 

1.33 

30.35,40 

260 

hitcnsity  distribution. 

1.27-1.3  Vi 

C  3 

276 

Efficiency  factor  for  extinction,  radiation 
pressure,  absorption  and  scattering. 

1.29(1-:U 

0-20 

278 

Variation  of  extinction  curves  if  the 
imaginarv  p-a-t  of  the  refractive  index  is 
varied. 

G.9-.69i 
8.18-1. 9fii 

00 

.7-1. 5 

283 

Extinction  curves  (showing  resonance 
peaTts) 

00 

3.41-1  94i 

0-3 

285 

Radar  cross  section  a  comjiuted  for 
backscattering  by  watei  drops  at 
X=3nim. 
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Some  representative  bistatic  cross  section  curves  for  a  perfectly  conducting 
sphere  are  included  here  ( Figures  2-4  through  2-20)  to  indicate  the  behavior  of  the 
sphere  as  a  scatterer.  These  were  computed  at  Air  Force  Cambridge  Research 
Laboratories  and  appear  in  King  and  Wu  |l95^  .  The  back  scattering  cross  section 
as  a  function  of  ka  lor  the  perfectly  conducting  sphere  is  also  given.  Figure  2-21. 
This  was  plotted  from  the  tables  of  Hey  et  al  |l95^  . 
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PERFECTLY  CONDUCTINO  SPHERE.-  kii  ^  3. 5 
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LOW  FREQUENCIES 

At  low  frequencies  the  scattering  function  1(9, 0)  can  be  'panded  in  series 
of  ascemiing  (positive)  powers  of  k  with  coefficients  which  are  functions  of  6  and  0. 
Associated  with  each  power  of  k  is  the  corresponding  power  of  a  parameter  Jt 
having  the  dimensions  oflei^ith,  and  since  f(d.0)  is  independent  of  r,  this  param^er 
must  be  a  characteristic  of  the  scattering  body  It  is  obvious  that  in  the  case  of  a 
sphere  the  parameter  is  the  radius.  For  sufriciectly  small  values  cf  \Jt  this  expan¬ 
sion  is  absolutely  convergent  (a  fuller  discussion  sf  the  convergence  properties  is 
given  in  section  3.4),  and  is  generally  referred  to  ar  the  Rayleigh  series  for  *he 
body  in  question. 

The  present  section  is  cittirely  concerned  with  this  expansion,  and  the 
purpose  is  not  only  to  dr*'  rmine  the  form  of  the  series  (i.  e.  the  powers  of  k  which 
it  conta:  but  also  the  precise  coefficients  of  the  various  powers 

In  section  3.1  the  scries  is  detained  directly  from  the  Mie  sol*ition  by  ex¬ 
panding  for  small  ar^piment  the  s{AericaI  Bessel  and  Hankel  functions  occurring  in 
the  solution.  la  so  doing  the  aim  was  to  set  down  explicitly  a  significant  nur  :ber  of 
terms  in  the  expansion,  and  presented  here  are  the  first  five  terms  in  the  exp.-!.'£ion 
fov  the  real  |»rt  of  1(9, 0),  together  with  the  first  four  terms  in  the  expansion  fur  : 

the  imaginary'  pert.  The  resulting  expression  for  the  scatterii^  junction  then  | 
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includes  terms  in  p’“.  where  p  =  ka 

In  section  2. 2  an  alternative  methofi  is  developed  wherel^'  the  Rayleigh 
series  is  obtained  directly  without  any  reference  to  the  Mie  soltitton.  aiai  without 
ever  having  to  solve  a  Ijoundary  value  problem  as  such  There  is  no  limit  to  the 
number  of  terms  which  can  be  calculated  in  this  way,  and  while  the  derivation  of 
the  higher  order  terms  can  become  tediosis.  ‘hs  labour  is  nc  worse  than  that 
involved  in  the  expansion  of  the  Mic  coefiTcieits.  In  addition,  the  calculation  is 
partially  self-checking. 

One  of  the  main  ath'antages  of  this  new  approach  is  the  promise  whirl:  »t 
holds  of  being  applicable  to  other  (and  more  general}  1/odies  for  which  the  exact 
Mie-type  solution  is  not  available,  but  even  with  a  spherically  stratiHed  sphere  it 
may  be  quicker  to  use  this  method  to  obtain  the  first  few  terms  in  the  Rayleigh 
series,  and  hi  section  3.3  the  leading  terrn  for  a  dielectric  coated  sphere 's 
calculsteu. 

3.1  Derivation  from  the  Mie  Series 

Since  tiie  exact  solution  for  the  sphere  is  known  in  the  for:;-  of  the  Mie 
senes  it  is  only  necessary  to  expand  the  radial  functions  for  small  o  :o  otHasa  the 
Rayleigh  scrics. 

The  coefficients  of  the  vector  wave  functions  for  a  iierfectly  cc  JveUr.^r 
spliere  are  given  in  eqns  (2-25),  and  using  the  fad  that 
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i  7  /,  7  2  91  4\  13 

2023  ^  y  ■  60  ^  "  ^O"  ^  J*  > 


j  0 


ol-i  20<ir:5>2 


1  s 


el4  <420P 


|5  P  {  1-  Tnjr  P 


153  2 

1540  ^ 


V  o(p^^> 


A  -  ^  M  r.  /  ^3. 

ol5  •  30(945P  ^  -  0(p  ) 


1%  i  _  11  g  l3« 

%ro  lirosT  > 

For  n?  and  B^j^are  CM  . 

The  above  sxpu.sions  are  sufficicEJ  to  specify  f  {0,  0)  correct  to  0{p\ 
but  rather  than  write  down  the  resultii^  series  for  arbitrary  6  and  0  we  shall 
concemr  ue  on  the  particular  cases  of  back  and  foiW'ard  scatterii^  (0  =  0  and  sr 
respectively),  hubstitution  into  eqns  (2-42)  and  ,^-43;  ihen  gives 


17  4 

900  *' 


S65I923  6  249170261  sl 

7938000  ^  "1875352500  ^  f 


1  -  6  J 

^ip  I  1+- 


I  i 

6804  "  I  ’ 
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f  (r) 


1_  3  J  1|3_  2  li83  4  670057 
2  ^  i  '  *  M  ^  ~  21W  ^  ~  793800 


2137  4 

&45(»*  ^ 


6  836^55737  3  I 

^  '  63762&25W  ^  J 

56639  6 

56700  ^ 

J 


C3-4) 


3.2  The  Alternativft  Meihod 

In  order  to  illustrate  the  ierm-l^'-term  technlqi»  for  deriving  the  low 
fi^uency  expansion  it  is  convenient  to  consider  once  agiun  the  imAlem  of  the  field 


teq;:S  {2-1)  )  incident  on  a  pcrfecUy  cCKiducling  si*crv. 


The  first  step  is  to  postulate  a  general  *.xr*--ssion  for  the  scattered  Celd 


and  the  trfivious  form  is  that  shown  in  eqn  (2-35).  Each  of  the  vector  wave 


i  f-inctions  M  ,  and  X  ,  involve  the  radial  distance  r  through  the  Hankel  function 
I  -oln  -eln 

I  aixi  its  derivatives  with  respect  to  fcr,  aiKi  conseqtnmliy  any  power  of  r  is 

I 

t 

*  always  aextompanied  by  a  like  tanver  ot  K.  Near  to  tte  surface  of  the  sj&ei^  r 

*  and  for  safOcieatiy  small  valitts  of  kr  (that  is,  for  stfiicisntly  hm  &eq^«;ies) 


h  (kr) 
a 


(kr) 


n+1 


(2n)! 

2“nl 


(3-:.) 


As  a  result,  all  the  components  of  N  ,  are  of  equal  order  in  the  i^ar  field  (in 

~eln 

contrast  to  their  behavior  in  the  fcr  Held),  while  ite  osmpouenls  of 

sa--ee  value  of  n  are  of  one  higher  order  If,  therefore,  the  prodwl  B  r  ,,  Is 

©II 

o  remain  finiie  in  ilm  scar  field  as  the  fz^iuency  c^cx^ases  irK^nnliely,  H  ts 
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necessary  that 


and  since  is  dimensionless,  k  must  be  associated  with  a  length  parameter 

3 

which  can  only  be  the  rad'"?  of  the  sphere.  Hence  B  ,,  --  O  (p  )  for  small  p,  and 

ell 

from  a  consideration  of  the  higher  powers  of  kr  in  the  expansion  of  jjcr  h^  {kr)J 
for  small  kr  we  are  led  to  write 

®ell  ^  <'^10  ^  ^  . ^ 

Similarly 

®el2  ^  ^  ^^20  ^^21  ^  ^22*^  ^23  . ^ 

and  so  on.  Any  of  the  above  coefficients  may,  of  course,  be  zero. 

For  the  prod'Jct  a  like  analysis  would  siiggest  that  the  expansion 

2 

for  A^^j  should  start  with  a  term  in  p  ,  but  by  choosing  instead  th.'  eypressio.a  for 

the  maioielic  field  neai  to  the  surface  (so  that  A  ,,  occurs  in  combination  with  N  „) 

oil  —on. 

2 

it  is  seen  that  the  coefficient  of  p  is  in  fact  zero  V/c  therefore  take 


^oll  ^  "l2''P  •^13''  ■  •  ' 

analogous  to  equation  (3-  6),  and  similarly 


(3-3) 


56 


THE  UNIVERSITY  OF  MICHIGAN 

3648-1-7 


At  the  surface  r  •  a  the  boundary  conditions  require  the  vanishing  o£  the 

tan iential  components  cf  the  total  electric  field,  and  substituting  the  expressions 
i  s 

for  E  and  E  ,  we  have 


X  + 


and 


f  /  \  8  f  \  r 

y  j  f A  ,  h  (p))  —  P  (cos5)+  (B  ,  — fph  (p)  I  )  ■“.  . 
t  \  oln  n  /  9b  n  ^  eln  n  n  Jj  sin  0 


■  0. 
(3-10) 


X  cos  ^  ^  ^ 


1/. 


i  A  ,  h  (o)i - ^ 

\  oln  n  /  cos  0 


\  P  (cos0) 
.»  n 


-  0, 


f  where 


X  =  e 


-ipcosb 


v-ipcosb)” 

n.' 


(3-11) 


(3-12) 


I  1  r  . 

{  Since  the  expansions  for  ’  (p)  and —  Iph  (p)  I  are  known,  the  coefficients  o..  and 
I  n  p  1^  n  J  xj 

B,.  in  t’  c  expansions  for  the  A  and  B  can  now  be  determined  hv  eqiintiag  to  zero 

il  on  en  .  .  & 

the  coefficients  of  each  penver  cf  o  in  equat'oas  (3-10)  and  (3-11). 

In  both  equations  the  lowest  power  cf  p  ic  p°  and  the  coefficient  is  made  up 
of  a  contribu'ion  from  the  single  vector  wave  function  and  the  static  teiiii  in 
the  incident  field  expansion.  To  this  oxrder  in  p  the  two  boundaxT’*  conditions  reduce 


to 


(I  V  i^j„)  -  0 
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j 

i 

and 

(1  cos  0  =  0 

giving 

<T  =  t 
*"10 

The  second  stage  in  che  analysis  involves  the  terms  in  fi.  Contilbutions 

from  the  two  furcher  wave  functions  M  ,  and  N  are  now  introduced,  together 

^oll  ”el2 

with  a  contribution  from  and  are  matched  to  the  second  term  in  the  incident 
field  expansion.  We  have 

-icos0  -  iOjgCose  +  +  ISig^gOos  0  =  0 

-i  cos‘°0  -  iOjQ  +  i^j^  cos  0  +  18  i^^^,  cos  20  =  0 

and  by  identifying  coefficients  of  like  trigonometrical  functions  in  each  equation,  it 
is  found  that 

^11  "  “10  '  T'  ^20  "  36 '  • 

Continuing  in  this  manner  the  variocs-  terms  in  the  expansions  of  the 
and  can  be  derived,  but  since  the  analysis  is  so  entirely  straight  *'orward  there 
is  little  point  in  including  further  stages.  Suffice  to  say  that  the  results  are  in 
accordance  with  those  given  in  section  3. 1. 

On  the  other  hand,  there  are  several  features  of  the  method  which  it  is 
desirable  to  point  out.  In  the  first  place  we  remark  that  the  analysis  at  each  stage 

! 
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is  to  some  extent  self-checkii.g,  in  that  the  nth  stage  provides  2n  +  1  s^if-consistent 
equations  from  which  to  calculate  2n-l  uiiknowns.  Moi'eover,  the  nth  stage  (which 
brings  in  contr:  jutions  from  the  wave  functions  j  — elu^  requires  that 


n-ll 


0 


'n-3  3  ’ 


and  a 


'n-2  y  V4  3'  “n-6  5 . 

all  have  the  value  zero,  aud  therefore  introduces  no  new  power  of  p  into  the 
expansions  for  the  corresponding  coefficients  of  the  vector  wave  functions.  In  fact, 
each  stage  yields  a  correction  term  to  tne  expar.sion  for  either  or  ir  <  n), 
but  not  both,  and  since  and  are  of  the  same  order  in  p,  two  successive 
stages  are  needed  to  give  a  new  order  of  correction  to  both  the:.j  coefficients. 

A  further  point  of  interest  concerns  the  real  or  imaginary  character  of  the 
a-j.  and  the  .  At  every  odd  stage  in  the  aiudysis  an  even  power  of  p  is  matched 
to  a  like  power  of  p  in  ti..;  incident  field  expansion,  and  from  eqn  (3-1?)  it  is 
stppaxe.tX  That  this  implies  the  matching  of  the  appropriate  a. .  and  0  .  to  a  real 

ij  ij 

coefficient.  In  contrast,  the  even  stages  produce  values  of  a,,  and  which  are 


pure  imaginary  and  hence 
real 


odd 


O',,  is  ,  .  if{i-j}  is 

ij  imaginary  '  even 


wiereas 


o  ,  real  .,  ,,  .  even 

3..  is  ,  ,  if  (i  -  j)  IS  , , 

ij  imaginary  odd 
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it  now  .‘ollcws  that  when  n  is  odd  all  even  powers  of  p  in  the  expansions  for  and 

iB  ,  have  real  coefficients  (odd  powers  having  pure  imaginary  coefficients),  and 
ein 

the  reverse  situation  holds  where  n  is  even.  This  immediately  determines  the 
power  of  p  in  the  expression  for  f(6,  0)  which  have  real  or  imsginarj'  coefficients, 
and  reference  toeq.i  (2-3?)shows  that  all  even  powers  must  have  imaginar5' 
coefficients,  while  the  odd  powers  have  real  coefficients.  The  fact  that  the  first 

Q 

imaginaiy  coefficient  is  0  (p  )  is  a  consequence  of  the  vanishing  of  o^  and  |3^, 

4 

v.’hich  thereby  removes  the  p  powers.  These  conclusions  are  confirmed  by 
eqns  (3-3)  and  (3-4). 

Our  final  remarks  concern  the  initial  stages  in  the  analysis  At  the  first 

stage  the  coefficients  of  p°  are  matched  and  this  requires  that  the  incident  field 

•‘ikz  i  i 

factor  e  be  replaced  by  unity,  so  that  E  and  H  are  independent  of  one  <  lother 

to  this  ai^roximation.  Moreover,  onlv  conlrihiitAo  a  of  order  p°,  and 

conseque  Uiy  this  first  approximation  has  produced  a  near-field  boundary-value 

problem  in  which  the  electric  and  magnetic  fields  are  nccoupled.  Although  the 

coupling  is  re-introduced  at  the  second  stage,  it  may  be  of  interest  to  consider 

why  the  initial  deco  upling  does  not  affect  the  derivation  of  a  complete  solution. 

The  first  stage  essentially  reduces  the  problem  to  a  static  one  and  gives 

I  only  (IjQ  =  i,  which  corresponds  to  a  simple  electric  dipole.  Tlius.  Iht  fivst  tiag- 

ignorcs  the  magnetic  dipole  contribution  to  the  scattered  field,  which  contribution 

is  of  the  same  order  as  the  electric  one,  and  to  obtain  the  electric  field  due  to 

i 

Uini  — .  . . .  6C  I.  I.  I.  I— ...  ■  —I—  ■  ■  ■  ■■ 
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the  magnetic  dipole  either  of  two  methods  can  be  adopted.  The  first  of  these  would 

require  the  corresponding  first  stage  in  the  solution  of  the  magnetic  field  problem 

and  the  subsequent  use  of  the  field  relations  to  determine  the  contribution  to  E  . 

In  practice,  however,  this  is  not  necessary  in  that  the  second  stage  in  the  solution 

of  the  electric  field  pi’oblem  re-introduces  the  coupling  between  the  electric  and 

magnetic  fields  and  brings  in  the  magnetic  dipole  contribution.  Two  stages  are 

therefore  necessary  to  complete  the  first  approximation  to  the  scattered  field,  and 

the  fact  that  no  magnetic  field  problem  as  such  has  to  be  considered  is  a  direct 

consequence  of  the  symmetry  between  the  expressions  for  the  scattered  61601.1*10 

and  magnetic  fields  in  terms  of  the  wave  functions  M  and  N 

-e  nm  -e  mn 

o  c 

I  On  the  other  hand,  if  only  the  first  term  in  the  expansion  for  f(9, 0)  is 

I 

I  required,  it  may  be  more  conveoient  to  replace  the  second  stage  by  the  first  stage 

« 

I  of  the  .correspondiog  magnetic  dipole  analysis,  since  this  may  prove  to  be  a  some¬ 
what  easier  calculation  (particularly  for  bodies  other  than  the  simple  homogeneous 
si^re).  In  this  case,  the  whole  analysis  can  be  oApicSesed  more  concisely. 

Taking  first  the  electric  dipole  prc^lem.  t!m  first  stage  is  to  match 
®ell  ^11  ’  cctor  ^  at  the  surface  r  =  a  using  the  boundary  condition 

n-CE'  +  E*)  -0. 

aru  since 

$  =  sin  d  cos  0  ^  -r  cos  6  cos  0  &  -  sin  00  , 
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consideration  of  the  tangential  components  of  shoH'S  that 


®ell  ’  TTTTTV"  for  small  p. 


— [pr,(p)] 


For  the  magnetic  dipole  problem  the  corresponding  stage  is  to  match  i 
to  the  unit  vector  y  using  the  boundary*  condition 

n  •  (H*  +  H®)  =  0. 
and  in  like  maiii‘.er  this  g'ves 


oU  '  2 


1  3 

for  small  p  . 


This  completes  the  analysis  for  the  two  near-static  problems  The  electric 
dipole  rnakcs  a  direct  contributicn  to  the  scattere'*  electric  Held,  and  according  to 
the  first -of  eqns  (2-35)  we  iiavc 


c®  i  E  N 


o  -ell 


(3-13) 


Similarly,  the  magnetic  dipole  contributes  direciiy  to  me  scattered  magnetic  field 
and  from  the  second  of  eqns  (2-35) 


J_  ::  V 

“  Xfi  *1 


c‘-^ll  ’ 


(3-14) 


''rom  which  the  electric  field  contribution  can  be  found  by  using  Maxwell's 
equations.  The  importance  of  this  derivatiou  lies  ir  the  fact  that  both  (3-13)  and 
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(3-14)  can  be  obtained  by  a^^aling  only  to  statics. 

!n  practice  however,  the  last  step  (use  of  Maxwet.’s  equations)  can  be 
avoided  by  substituting  the  expressions  for  A^jj  and  directly  into  the  first  of 
eqn  (2-35).  We  then  have 


which  represenlb  the  combined  contribution  due  to  the  electric  and  magnetic  dipoles, 
and  the  correspondii^  far  field  expansion  is 
_ikr 


O  XT 


3^  ^  9)  COS  0  5  -  (1  s  —  COS  d)  sin  0  ^  * .  (3-15) 


j  3.3  A  Dielectric-Coated  Sphere 
( 

I  As  an  example  of  how  the  above  method  is  used  in  a  non-trivial  problem,  we 

» 

shall  here  derive  the  leading  term  in  the  Bayleigh  solution  for  a  coated  sphere 
(see  seen*  n  2.2). 

Consider  a  perfectly  conducting  sphere  of  .  a  which  is  covered  with  a 
Ia3rer  of  dielectric  of  thickiKss  d.  The  permittivity  and  parmeability  of  tte 
dielectric  are  f  and  ^  respectively;  the  conductivity,  however,  is  zero  totherwise 
the  Rayleigh  solution  is  i‘te  same,  to  the  first  term,  as  for  a  perfectly  conducti^ 
sphere  of  nulius  a  +  d).  The  whole  is  immersed  in  a  homc^neous  iso^^cpic  mediuit: 
which,  for  simplicity,  will  be  regarded  as  free  space. 
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FIGURE  3-1 


To  determine  the  Rayleigh  solution  it  is  sufficient  to  solve  two  static 


problems  and  then  match  these  results  to  expressions  involving  spherical  vector 
wave  functions.  In  the  first  static  problem  the  incident  field  is  merely 

=  E  X  (3-J6) 

o 

and  ihe  task  is  to  obtain  the  scattered  (electrostatic)  field  which  this  excites.  Tljc 
second  problem  is  analogous  in  that  the  incident  field  is  here 


=  -  H  9  . 


<3-ir, 


so  i’*at  a  magn^ostatic  field  is  now  involved. 

The  most  general  solutions  of  Laplace's  elation  are  of  the  form 


_  w 

5  =■  r°  P™  (cos  m  0 

i  n  sin 


^  P  (cos  0)  .  m  p 
n  sin 
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and  if  ^  is  legarded  as  a  static  potential,  the  corresporniing  field  can  be  found  by 

taking  the  gradient.  In  the  region  outside  the  sf^re  the  scattered  field  must  be 

(2) 

expressible  in  terms  of  £  alone,  but  in  the  dielectric  coating  both  types  of 

^  e  mn 
o 

potential  will  occur. 

Let  us  take  first  the  electrostatic  problem  in  which  the  fr.i  ident  field  is  given 
by  eqn  (3-16).  Since 

X  =  V  J 

ell 

we  have  immediately  that 
i 

E  -  E  V  §  (3-18) 

®  en 

18  * 

If  E  is  the  total  electrostatic  field  in  £r^  space,  so  that  M  -  E  •  K  ,  and  if  E 
i 

[  similarly  d^iotes  the  field  in  t.*!e  layer,  the  boundary  conditions  at  the  dielectric 
interface  (r  =  at  ui  ai.- 

^  /N  E  =  A  E* 

and 

€  a  E  =  e,  n  ■  E 

o  -  1  - 

At  the  surface  {r=a)  of  the  perfectly  conducting  sphere  the  omy  (xmdHitm  is 
n>^E’  -  0  . 

I  !n  view  of  the  &  and  0  depei^nce  implicit  in  the  exjuression^qn  C3-]8^for 
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the  incident  field,  it  is  apparent  that  the  boundary  conditions  can  be  satisfied  by 
choosing  the  following  expressions  for  the  secondary  field: 


E®  -  A  V  I 


E*  =  B  VO  +  CVO 

eU  ell 


<3-20) 


where  A,  B  and  C  are  CfMistaitts  as  yet  undetermined.  The  boonditxy  cffisditioos  now 


E  +  - =  C  +  -  , 

®  (a+d)"*  (a+d)^ 


g  1  E  -  — A 
®  1  ®  (a+d)^ 


^  *  €■  C- 
i 


0. 


from  whiL-r  we  obtain 


A  ^  -  E  b' 
o 


J  b^2a^--^(b^-a^) 

I  b^  +  2  a^  +  2  ^  (b^  -  a^) 


<3-21) 


where  b  =  a^.  If  d  «  a  so  that  powers  of  d/a  higher  than  the  first  can  be 
neglected,  the  expressicm  for  A  be<mmes 
_  .  .  ..3  ^  d A 


A  .  («  —  f)- 


<3-22) 
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which  differs  from  the  result  for  &  perfectly  coitduciing  sphere  of  radius  avu  only 
in  the  presence  of  the  muUtplyiug  factor  A- 3 


Frcin  eqos  (3-19)  and  (3-221,  the  field  which  is  scattered  into  free  space  is 

e  (2) 

$ 

ell 

and  tt«  next  step  is  to  match  this  to  the  limit  ot  a  non-static  solution  at  low 
freqeeuctes .  Si  ncc 


E®  =  -  (a+d)^  ^1-3  5 


2  sin  e  cos  p  i-  cos  6  cos  0  ©-sin  0 


?!■ 


consioeratiai  of  the  vector  wave  functions  M  and  N  shows  that  for  X»  r, 

-emn  -e  am 

o  o 


V  I 


(2) 


i  ,, 


eu 


and  hence  m  tbt  near-static  litrit 


E®  =  !  E  (  1-3  ;  N 

-  o  \  e  a  /  -« 


*-^U 


(3-23) 


The  far  Jield  is  now  obtained  by  inserlins  the  firs*  terms  of  the  asymptotic  expan¬ 
sions  of  the  mdiai  H.tnkel  sur.ctions  for  large* kr,  srtd  this  gives. 


ikr 


E®'"E^  *kr”  *‘-"{3-fd)^l-3  —  (cos  9  cos  0  3-sin  0  0)  . 


S7 


THE  UNIVERSITY  OF  MICHIGAN 

3G48-1-T 


Eqn  (;j-24)  represents  the  k  contribution  to  the  far  field  arising  from  the 
electric  dipole,  but  it  is  not  the  only  contribution  of  this  order.  There  is  in  addi¬ 
tion  3  term  produced  by  the  magnetic  dipole  and  this  is  most  conveniently  (Stained 
by  considering  a  second  static  problem. 

Tile  incident  (magnetostatic)  Held  is  no«- 

y'-V  • 


uhich  can  be  writtmi  as 


h'=-h 

—  o 


(1) 

oU  • 


(?-25) 


and  the  task  is  to  find  the  scattered  field  subject  to  the  boundary  conditioi.s 

ji  n-  H-fcn>  H 
o  ~  ~ 

n*B:  n^'H 


at  r  =  a+d,  and 


i>  •  ts  =  o 


at  r-a.  The  fc-m  of  emt<3-25)  leads  us  to  tulc^t  the  following  expressions  for  the 


1 

I  where  A,  B.  and  C  are  cmistants  as  yet  undetermined,  aixi  the  boundary  conditions 
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then  gi'.-e 


v{-V 


A  .  B 

-H  4- - --r-  ^ 

o  (a+d)^  ^  (a+d)^ 


C-^-0 


Solving  for  A,  w  have 


f 


H  b-^  [  (b^-a-") 

T - 2_  i - 2 - 

“  2  I  ^3  3  ti-i  3. 

I  A  +a  +y  (b  -a  ) 

»■  It 

abere  b=a+d,  and  if  d<<  a. 


(.?-2S) 


i)- 


(3-27) 


This  differs  from  the  result  for  a  perfectly  eoMlacUr^  sphere  of  rsKiius  a+d  only 
}r«  ifee  presene-  -«  nsuUi^ying  factor  ^1-3  t)  • 


'^ne  field  which  is  soitered  in!c  free  space  is  mnv 

t  y-  S', 


M 


t(2) 


2  , 
ai hy  maicht:^  to  the  vector  uavc  functions  for  A  »  r.  we  have  in  the  n  tir- 
stiiic  Umii 

ii 


n  ®=  i  Y  (l-3  ^  I)  . 
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The  correspondinp;  electric  vector  is,  by  using  Maxwell’s  equations 


S  *^0  ,  3,  ,  ,.3  / ,  ..  M  d  \ 
E  =  -  Y  k  (a+d)  1^1-3 


%11 


(3-29) 


and  in  the  far  field  this  becomes 

ikr  ,  3,  ^  ,,3 
..s  „  e  k  (a+d) 

E  "E  —  — - - 

o  kr  2 


{  1-3  ^  f)  (cos0  0-c  >f  6  COS0  0).  (3-30) 


The  complete  first  term  in  the  Rayleigh  expansion  for  the  scattered 

I 

’  electric  field  is  obtained  by  adding  the  ''ontriuuiions  represented  by  eqns  t3- 
and  (3-30).  Tlic  required  solution  is  therefore 

^)cos0]0in0^|  . 


(3-31) 


from  which  the  scattering  function  ca  t  be  determined  if  so  desired. 

3. 4  Convenience 

‘  -  previously  remarked,  the  Rayleigh  serie<?  is  a  convergent  representation 
for  sufficiently  small  values  of  k  /  ,  and  in  ■'  y  application  of  the  above  rc.sults  the 
actual  ralius  of  convergence  is  then  a  matter  of  some  importance. 

To  see  how  the  convergence  arises,  let  us  write  the  .scattered  electric 
field  in  the  for.m 


a  iki)  f  (r,O,0) 
n  ~n 


«3-32) 


■10 
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where  f  (r,  9,0)  is  a  vector  function  of  the  coordinates.  The  series  <  n  Ihi.  right 
~n 

hand  side  is  absolutely  convergent  for  all  values  of  k  ^ ,  the  function^  f  {r,  6, 0) 

~n 

being  bounded  as  functions  of  n.  Each  a  (ki)  can  be  expanded  in  a  series  of  posi- 

n 

tive  povers  of  kJ  in  a  neighborhood  of  the  origin  of  *he  complex  k£  plane,  and  is 

therefore  an  analyiic  funedon  of  k/  within  this  region,  by  rearianging  the  terms 

in  eqn  (3-32)  we  then  have  a  representation  for  as,  an  expansion  ..  powers  of 

k£ ,  which  expansion  converges  within  the  least  circle  of  convergence  of  the 

individual  a  . 

n 

If  the  functions  a  are  now  identified  wiT”  th<-  coefficients  A  ,  and  B  , 
n  oln  c!n 

I  in  the  vector  wave  function  expansion  for  a  perfectly  conaucting  sphere,  it  is  a 
simple  task  to  determine  the  appropriate  radius  of  convergence.  From  eqn  (2-25) 
it  is  apparent  that  the  only  singularitits  of  the  A^^^  and  are  poles  at  the 

I 

I  zeros  of  the  spherical  I.'aiikel  function  or  its  derivaiue,  and  the  location  of  l.iese 
;  zeros  is  such  that  the  singularity  nearest  to  the  origin  is  providaa  by  one  of  the 
smaller  ...^aes  ot  n.  For  n  =  1  we  have 

hi(p)=-^  (1  +  ^ 


showing  that  A^jj  has  a  pole  at  p  =  -i  and  has  poles  at  p=  -  ^(1±  i,*^! 
Accordingly,  both  A^jj  and  are  infinite  on  the  urit  circle  and  since  all 
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the  higher  coefficients  are  regular  inside,  the  entire  Rayleigh  series  must  con¬ 
verge  for  |p|<l.  The  fact  that  a  singularity  exists  for  which  |p|  =  1  implies 
that  tie  series  does  not  converge  outside  this  region,  and  consequently  the  Raylei^ 
series  for  the  perfectly  conducting  sf^ere  converges  only  for 

ka<l.  (3-33) 

From  the  above  discussion  it  is  obvious  that  the  convergence  is  determined 

solely  by  the  coefficients  A  ,  ajid  B  ,  ,  and  any  change  ’a  these  by,  for  example, 

oin  ein 

a  modification  cf  the  boundary  condition  may  affect  the  overall  convergence  of  the 
Rayioigh  series.  To  illustrate  this  point,  Ic,  us  consider  the  case  in  which  the 
boundary'  condition 

E-(a.  E)n  =  q  Zn^H  (3-34) 

is  imposed  at  the  surface  of  tlie  sphere.  Here  n  is  a  unit  vector  normal  drawn 
outwards  ft  om  ♦he  sphere,  q  is  tlie  reciprocal  of  the  complex  refractive  index  of 
the  I'-aierial  of  the  sphere  relative  to  free  space,  and  Z  is  ♦iie  intrinsic  impedance 
of  free  space. 

Eqn  (3-34)  is  the  usual  impedance  boundary  condition  and  is  only  accurn^c 
to  the  first  order  in  q  .  The  physical  situation  therefore  requires  that  r  he  assumed 
small  (q=0  for  infinite  conductivity),  though  there  are  circumstances  unde.*  which  a 
I  physical  significance  can  be  attached  to  eqn  (3-34)  even  when  q  I®  not  email  Ci...  - 
pared  with  unity.  On  the  othei  hand,  if  the  problem  is  merely  regarded  as  a 

.  ~  I  . .  II  I  .......r,  ■  n....—  72  ■  -  - 
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mathematical  one,  it  is  a  tri/ial  matter  to  show  that  the  boundary  condition  is 
precisely  satisfied  by  a  scatt;!red  field  of  the  form  (eqn  (2-2  6)  )  with  coetficients 


where 


A  =  fl  (p.  irj), 
oln  f* 


'  V”- 


u  (p,iy)= 

n 


-(-i) 


n  '2n-i-l 
n(n+l) 


h^(pKwi  [ph__(p^' 


(335) 


(3-36) 


For  a  fixed  value  of  V  ,  f2  (o,  iV)  is  a  function  of  p  an-  can  be  ext)anded 

n 

in  a  convei'gent  power  series  within  some  neifdiborhood  of  the  origUi  p=0.  The 
circle  of  convergence  depends  on  n  and  as  it»  the  cr  e  of  a  perfectly  conducting 

{ 

»  sphere  the  least  circle  is  pretiJed  by  A  „  and  B  ,,,  that  is,  by  The 

j  oil  ell 

denominator  in  the  e.vpression  for  f2i(p,  iV)  Is 


J 

•  which  vanishes  when 


These  are  two  genuine  /.eros  except  «4ien  V  =0,  in  which  case  the  second  zero  must 
he  •liscounted.  If  this  case  is,  for  the  nicment,  excluded,  i‘  follows  that  the 
Rayleigh  expansion  for  E  converges  only  for 
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When  Visa  general  complex  constant  an  explicit  form  for  |  P}  |  or  |  p2  ( is 
ditficuli  to  write  <.iown.  ll!  V  j»l  or  j  V  |«1,  however,  me  equations  for  Pi  and 
P;  simplify  considerably,  leading  to  a  more  compact  statemcn  of  the  convergence 
region.  Thus.  for|V  j»l. 

and  for  I V  |«  1, 


Pi — i(l+V) 

Pa- sV. 

Accordingly,  for  sm.iil  q  ihc  convcigcncc  legion  is  specified  by  the  zero  p-j  for 
the  coeffinien*  A^jj « .  replaced  by  q  )  and  is 

ka<  lql«l; 

similarly,  for  la.  ge  q  the  convergence  is  determined  by  the  zero  p?  for  the 
coefficient  B  ^  ^  <  V  replaced  by  l/rj)  and  is 


In  both  cases  the  region  of  convergence  is  appreciably  reduced  in  comparison  wi‘h 
that  for  a  perfectly  conducting  sphere,  and  can  become  infinitesimally  small.  We 
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observe,  however,  that  for  r,  =  0  the  zero  which  is  dictating  the  convergence 
disappears.  Th  i  zero  p  i  then  becomes  important  and  leads  to  the  result  given 
in  eQn(3-33). 

If  17  is  neither  lu2  go  nor  small  conipared  with  unity,  the  boundary  condition 
is  of  doubtful  validity,  but  it  is  still  of  interest  to  examine  the  convergence  of  liie 
Rayleigh  expansion  when  17-0(1).  Both  Pi  and  P2-VOO  asV-*-!  and,  indeed,  for 
V  =  17  -  1,  Pi  and  P2  are  inlinite.  'I  his  can  be  confirmed  by  looking  at  the  e.xpies- 
sion  for  fiifp,  i).  In  this  particular  circumstance,  however,  CyCp,  i)  imposes  a 
finite  radius  of  convergence  which  now  becomes  the  »ir.portant  one,  and  from  an 

I  examination  of  n2^P>  Oiat  the  Rayleigh  expansion  converges  only  for  ka<.3. 

On  the  other  hand,  note  that  17=1  corresponds  to  a  sphere  whose  impedance  is  that 
of  free  space,  and  this  is  certainly  a  body  for  which  the  impedance  boonUary  con- 

S 

I  dition  may  he  e.xpected  to  fail.  Nevertheless,  the  result  does  suggest  that  if  She 

i 

exact  boun-iary  conditions  were  used,  the  radius  of  convergence  may  be  greater 
than  unity  in  the  case  of  a  very  diffuse  sphere,  an'’  -  -.tu  !'  of  the  coefficients 
A  ,  and  B  in  Stratton  |l041,p.  3(»|  gives  additional  confirmation  of  this. 

Returning  now  to  the  previous  exanipie  in  which  17  is  large  or  smull.  tlie 
fact  that  a  marked  reduci.on  in  the  radius  of  convergence  of  the  Rayleigh  e.xpansion 
c' ccompanies  the  introduction  of  even  a  slight  impedance  into  the  sphere  is, 
physically,  rather  surprising,  and  suggests  that  Uic  usefulness  of  the  Rayleigh 
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approximation  is  limited  to  perfectly  conducting,  or  hi^ly  transparent,  bodies. 

The  discontinuous  change  in  convergence  between  the  cases  q=Q.  and  r]  ^  0  is 

dut.  to  the  fact  that  there  is  no  e.xpansion  for  A  ,  or  B  .  which  is  uniform  in 

oln  clr 

n-  tssentially  each  coefficient  involves  a  factor  of  the  form  p/ p-q,  and  for 
i7  ^  0  this  can  only  be  e.xpandcd  in  a  scries  of  positive  powers  of  p  when 
1  P  1^1  dl  •  Accordinglj'.  any  attempt  to  approximate  the  expressions  for  the 

of  (K/T)  will  be  limited  by  this  same  condition, 
somew’hat  disguise  the  fact. 


coefficients  by  neglecting  terms 
even  thou^  the  final  result  may 
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THE  WATSON  TRANSFORM  AND  CREEPING  WAVES 

The  problem  of  the  diffraction  of  electromagnetic  energy  by  a  perfectly 
reflecting  ^here  for  which  ka  is  sizeable  was  made  tractable  by  Watson  |^18> 

I919J .  Watson  found  a  transformation  of  the  Mie  series  -  the  Watstm  transform  - 
which  resulted  in  a  mudi  mo.e  rapidly  c<mvergent  r^resentatimi  of  the  solution. 
Mucli  later  Fock  ^1945, 1946j  and  Franz  ^35^  initiated  a  further  analysis  and  gen¬ 
eralization  uhlch  inr'icated  that  the  functional  form  of  the  Watson  solution  was 
ai^licable  to  problems  involving  other  convex  snapes.  The  mathematical  counter¬ 
part  of  the  extensions  of  Fock  and  Franz  is  found  in  the  work  of  Langer  ^193^  and  his 
followers  which  was  essentially  completed  for  this  application  in  the  1930's.  The 
more  general  approach  has  led  to  the  presentation  of  the  results  in  terms  of  certain 
"universal  fupcHons"  which  have  been  extensively  computed  and  tabulated  under  the 
direction  ji  N.  A.  Logan  |l950  ,  In  our  development  we  follow  the  appros.ch  of  Logan 
and  his  co-workers  j_188lj  . 

4. 1  The  Field  on  the  Surface 

We  now  compute  the  magnetic  lields  induced  on  Un;  suriace  of  a  perfectly 
ccmducting  sphere  by  pla;ie  electromagnetic  waves.  If  Uie  incident  electric  field  is 


Uiven  by 


-i  „  ^  “Ikz 
E  =E  X  e 
—  o 

-1  ..  ^  “>kz 

h  =-H^ye 
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the  magnetic  field  on  the  surface  is 


H  sin 

-Ir 


dP~^  (cos  e) 
n 


t"W) 


P  ^  (cos  0) 
1  n 

)  sin  e 

n 


H^cos  0  ro  1  ^P"  (cos  6) 


(ka> 

n 


(cos  0)1 


n 


where  wc  have  made  use  of  tiie  results  of  Section  n,  the  relati'^:; 


P  (cos  fl)  =  -  — - -  P  ^  (cos  6)  (4 

°  n(ii  +  i)  ** 

(1  2)  (1  2) 

and  the  notation  5  ’  (x)  =x  h  (x),  ^  (x)  =x  j  (x).  Foi-  convenience  «e  also 


n  n 


ma\e  the  substitution 


6— ►a  =  x-9 
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which  results  in 


sin  B 


sin  a 


(4-4) 


3P  ^{cos  0)  dP  "(cos  O') 

-5 - =  (-)“  “ 


36 


Using  equaticn  (4-4),  equaU<»  (4-2)  becomes 


H^sin  ^ 


0  ka 


.  on 


in»y2 


(2IH-1)  e  '  (— y 


a*l 


JT^^^'da) 


dP  (cos  a) 

D 

Jo 


5^^^'(ka) 

n 


P  ^(cosor) 
n 


sin  a 


(4-5) 


H  cos  f 
o 

u  -  I - 

p  ka 


<D 


Z-in*^^2  nl  i 

(2irt-l)  e  (-)“  I 

n«l  I 


3P  ^(rostt) 
n 


L  S 


S  (ks) 


i 


-n 


P  (cos  a) 
n 


sm  a 
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We  nov  rewrite  the  sums  over  the  integers  in  cq^.  (4-5)  as  sums  over  the 

odd  half-integers,  letting  a  =  7/-^  , 

H.sin0  T—  -i(2/-V2)'/2 

2Ve  (— ) 


xa 


i 


a--  f"v2"“’ 


4\«“' 


ll\^^'cosa) 

sina 


(4-6) 


«„•  ‘ 


H  cos  0 
o _ f 

ka 


% 

z. 

5-=  3/2 


-i{V-yz)^l2  V-\ 

2Ve  {-i 


1 


(I)' 


dP,.  1/  (cos  a) 
^-^2  + 


P  \  (cos  or) 

JL-^12^ 


(ka) 

^-72 


sin  a 


For  la*er  use  we  note  that  since 

p'  (x) =  0 
o 

the  te^'rrs  for  v  =  ^f2  could  have  been  included  in  these  sums. 

The  summands  in  eqn.  (4-6)  are  regular  functions  of  V  in  a  strip  along 
real  a.\is  so  that  we  can  write  the  sum  as  a  coitour  integral  about  the  positive  real 
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T  K-i  C»«nniac<fl 


jsinco  fOj  vr  has  siinplv  jioiss  at  the  odd  jialf- integers.  S^iecilically  the  eqns.  (4-6) 


I  become 


IF  sin  0  f  i/dp  -ii 

j,  ,  j  _2 - 1  - g 

if  ^2-  I  COS  P7 


•i{p-k)^h 


e  Ip) 


4  V 

j 

I  We  have  written  U»e  terms  in  square  bracScets  in  eqn.  (•!-{?)  as  B  and  ^ . 

I  IV^^2 

t  W'e  examine  the  lem.s  t  BiP}  and  acte  that  these  are  even  functimis 

I 

of  V.  Sera-js*-  iF.e  rsma.rakr  o!  the  integrand  is  odd  the  integrand  is  sa  o*^  ianction 
of  if.  WT  mslder  the  contour  C  in  Figure  4-1  snd  note  that  the  lower  path  givnis 

-plane 


HGURE  4-1:  THE  CONTOUR  C 
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an  integral  of  the  form 


oo-ic 


■‘7 


dvO{V) 


{4-9) 


0-l€ 


ttlicrc  VIC  wriU:  Ofv)  for  the  odd  functions  of  v  .  If  ac  reflect  the  conUMr  m  toe 
origin  a-e  find 


I  = 
e 


-a>Hc 

/ 

Oric 


dH  0(1') 


(4-10) 


However,  interchanging  the  limits 

(H-k 


dV0(l>) 


(4-U) 


-oo-rjc 

and  add.  ig  tl:e  coatributior.  f''  ♦he  upper  path  the  t«^  integral  is 

oprie 

dVOiV)  H-U) 

/ 

-ojf-k 

In  the  sequel  ae  will  sun>ress  the  k  in  tim  limit  with  the  onderstandiiK  that  tiie  rath 
is  to  run  just  above  the  real  axis. 

To  evaluate  the  int^trals  of  the  fmmi  of  eqn.  (4-12>  we  need  'X>  examine  the 
r-tegrand  in  smne  detaU.  Hie  Qrst  thiag  aote  i‘*  tii.it  tee  o)u.^  «.iids  have  simple 
poles  at  the  zerra  of  (la)  aad  (ks)  Aed  loese  lie  ia  the  first  qoadrant. 

>"/2  Jf*V2 
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This  analytic  behavior  s^i^ests  that  if  the  contour  (-a>,  co}  can  be  closed  by  a  semi 
circle  in  the  u;^r  half-plane  s^'e  can  evaluate  the  inle^.r3l  in  terms  ol  tnese  poles. 
This  is  indeed  the  case  umier  certain  circumstances  as  'ae  will  show  below.  In 
contrary  case  we  will  evaluaU  tnc  integral  by  the  method  of  stationary  phase. 

.4s  a  preliminary-  to  our  examinati<»'.  of  the  integrand  we  defire,  after  Logan 


fl96ll  ,  the  fiinctimss  9)  end  "24  S)  by 

L  J  m  3- 


<— >“  2P  “  {cos  fi)  5  {P.  £>)  +  9) 

V-/-)  m  m 

ttiiere  thesa  funcxicHis  have  the  asymptotic  beha%ior 


- ^  J-T^  9) 

m  Ij  sin  9  m 


(4-13) 


C4-14) 


for  |2;|-^aD  and  0<  e<  s  .  ExpUcitly 


ji).  [X:  i 

'^zsinS® 


2 


(4-15) 


1  1  ^ 

F  (-+  m.  T  -m,  3^+  i; - ) 

-  2  siin  fi 
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r* V2)  -  ,  -  r/  . 

*  2  '-i(z/0-  I  -m  72  ) 

e'7i/0)=~7=; - — - e 

m  P(-m-l)  \  ,r  sin  0 


.  -i0 

11  \ 

.  F,  (-T  m,  -  -  m;  V+  1;  - ) 

2  1 ' 2  2  2  siii  d 


(4-15) 

cont. 


We  have  remarked  that  the  contour  nms  above  the  real  axis  so  that 

Im  1/  >  0; 

hence,  we  can  make  the  convergent  expansion 

„  Ii/TT  ^  /  2jr  iVl 

sec  VTt  =  2e  y  (— )  e 

i  =  0 


(4-lC) 


i  i  the  form  of  E  ’  in  eqn.  (4-15)  we  make  the  following  observations: 
m 

(2)  -ii/ir  +1(^2  + ns*’/  Mi 

E\:^2^,^)=e  ‘^’"e  ^  e'V.  sr-a),  (4-17) 

•n  m 


(1)  (1) 

e  e'  '(•J-'.a)  =E  {.V.  a+J’jr). 
m  m 


(4-18) 


where  eqn.  (4-18)  is  derived  from  the  fact  that  the  hipergeometric  furotions  are 

periodic  in  c-  with  the  period  ir  so  that  the  continuation  of  the  E^^^  in  fne  ^-variable 

m 

iio  "  25r  is  determined  by  tlic  exponential  alone,  provided  we  take  the  radical 

I  (sin  v)  to  be  Isin  i/l"  '  .  Putting  these  results  together 

t  c  ...  ...  ■» 


•  seeVTT  P  (coso) 


2x  (i+D-  0)  -iEj  V.2jr.('+0)j 


(4-19) 
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where  we  have  returned  to  the  variable  8  =v-a. 

The  opei-ations  d or  1/sIn  6  on  P  \  do  not  essentially  affect  the 

behavior  of  P  , ,  as  a  function  of  for  1,  Im  V  >  0.  Hence  we  have  from  the 
Ty-Vo 

/I) 

asymptotic  form  of  the  E|‘  that  the  dominant  term  of  e<iu.(4-19)  is  of  the  form 
e  so  that  the  dominant  term  of  the  integrand  will  be  e  Therefo*’e,  for 

ImV  >  0  the  integrand  will  be  a  decreasing  exponential  in  provided  8  >  t/2. 

For  this  case  the  contour  can  be  closed,  the  semicircle  contribution  vanishes  ,  ;ind 
the  integral  is  given  by  the  residues  of  the  zeros  of  and  provided  the 

rest  of  the  integrand  remains  bounded.  This  is  indeed  the  case  except  on  the  locus 

of  the  zeros  of  (ka)  and  ^  (ka).  It  can  be  shuwn,  however,  that  a  path  can 

'"•''2  0)  -I). 

» be  found  between  any  two  zeros  on  which  the  functions  1/?  and  1/5'  '  remain 

: 


I  bounded. 


The  U;i  asb  otlier  Utau  uie  uuininant  one  in  eqn.  (4-19)  satisfy  the  convergence 
condition  on  the  scmicircie  for  all  values  of  8.  It  is  this  behavior  under  the  decom¬ 
position  (4-19)j  that  was  the  basis  for  the  "creeping  wave  •  analysis  of  Franz 


[l95^  . 


We  consider  the  behavior  of  an  integral  of  the  form 


GO 

1=  J  t>iV  e 
-CD 


m 


(4-20) 
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where  we  write  £  for  either  ?  or  .  Now  from  the  remarks  above  we 

//-I/2  vJk 

have  for  jr/2  th4it 

-it/  V2 

00  n  (n 

e  E,  {2/,W 

n  1 


I  =  2jri 


z 


(4-21) 


9V 


where  f(2/^)  =  0.  Since  the  first  zeros  of  f,  in  either  case,  occur  for }  V|  ^  ka  for 
ka  large  we  can  use  tiie  asymptotic  forms 


?'  ^  (ka)  =m  |wj(t)  -  ^  m  [(t^+9  )Wj(t)-4twj(t)J  +...J  (4-22) 


\s4iere  we  write 


1/3 


m  =  (ka/2)‘ 


t  =  -  (i/-ka) 
m 


(4-23) 


and  w^(t)  is  the  Airy  fimcUtm 

Wj(t)  =  {IT  (Bi  (x)  +  IAl  (x)  )  (4-24) 

Now  we  make  a  further  approximatitm  in  the  integrand  of  etpu  (4-20),  We 
give  all  slowly  varying  fimctiona  of  tf  their  values  at  2^  =  ka  and  remove  them  from 
the  integral.  This  we  do  since  for  t(/^it/2  ,  and  ka  and  hence  |  t)j^|large  encu^ 
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the  residue  aeries  [eqn.  (4-21)  ]  converges  with  sufficient  rapidity.  Again  uaing  the 
asymptotic  form  for  under  the  condition  \t^\  sin  tp  1,  we  have  in  this 

approximation  two  Integrals  which  we  write  as 

i?t 


«€)=  -pr 

iir 


,  CO  i5t 

=  f  - — 

It  j  Wj(t) 

J  -CO 


dt 


00 


(4-25) 


g(?)  = 


1  f  , 

]  w’d) 
•'-CD  ^ 


j  whcn  e  we  put 


€=m(  0-^). 


We  now  approximate  the  fields  on  the  surface  for  Q  >  jr/2  ,  ka  sin  0  »  1 


fcv 


00 


H. ‘  -  J  4 
d  o 


/=0  ^ 


(4-26) 


1=0  ^ 


where 


»  ni(2jr  /  +  fl--^)»m0^ 


5j|  =  m(2Tj^  +  6)  = 
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and,  as  above, 


m  =  {lui/2)^^^. 


The  next  region  on  the  surface  we  coasidei-  is  6  <  t/2.  Here  we  need  to 
further  decompose  ihe  region  since  the 'transition  from  0  o  to  the  shadow  boundary 
B  Ttl2  ,  is  accomplished  by  means  of  two  different  represent,  tjons.  First  near 
0  =  0  we  evaluate  the  fields  given  by  the  Integral  representations  of  the  form  of  eqn. 
(4-20)  by  a  saddle  point  method  for  the  first  term  in  the  expansion  ^eqn.  (4-19)^^  . 

The  result  for  this  leading  term  is 


-H  2  cosO  sin  0e 
0  o 


-ika  cos  6 


i  sin^  0 


5  8in^9-  8ln^9 


2ka  co8^9  2(ka)^coB®  6 


+  . . . 


(4-27) 


„  o  rt  -ika  cos  fl  1 
H0  -  -H  2  cos  0  e 


1  - 


1  sln^  Q 


9  sin^e  -  sln^O 


2ka  oott^O  2(l:a)^  cos^  6 


+  . . 


Here  we  "om;  rk  that  the  leading  terms  in  eqn.  (4-27)  are  Just  the  geometric  optics 
fields. 

To  bridge  the  gap  between  the  shadow  boundary  and  the  optics  fields  In  eqn. 
(4-27)  we  note  that  asymptotically  as  -oo  the  functions  f(C)  and  g(C)  in  eqn. 
(4-25)  go  to  the  correct  leading  term  of  the  optics  field  in  eqn.  (4-27).  iVe  can  then 
MSP  eqn.  (4-26)  provided  we  make  a  substitution  In  the  arg’oment  so  th^t  we  get  the 


correct  phase.  This  la  simply  to  let 


m(?  -  ^  »  m  sin  (9  -  | ) 


(1-28) 
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in  the  loading  terms  of  eqn.  (4-26)  for  Q  ^-n  and  Here  the  choice  of 

0^  is  somewhat  arbitrary.  To  make  the  choice  of  specific  is  not  meaningful  so 
we  specify  the  range  60°J  so  that  the  regions  of  eqns.  (4-27)  and  (4-28) 

overlap. 

In  elthei'  of  the  cases  we  also  remark  that  the  terms  In  eqn.  (4-26)  for 
X  1<  2,  ....  wiJl  also  appear  Just  as  before.  However  the  terms  near  the  caustic 
In  the  lit  region  diverge  as  l/^sin  0  as  0-»O.  We  will  find  a  bounded  represen¬ 
tation  for  these  when  we  treat  with  the  same  behavior  at  the  caustic  In  the  shadow 
region,  0  ■  y . 


2  Wo  will  now  find  a  representation  for  the  fields  in  the  region  near  the  caustic 

In  the  shadow  where  ka(r-0)  Is  small.  After  Logan  |196lj  we  make  the  physical 

argunjent  that  the  terms  in  eqns.  (4-26)  of  the  form  g(5^  )e^**®^/jslnT  describe 

v/aves  which  diverge  from  0  =  jr  while  terms  of  the  form  g(5^  )e  l\  sin  0  describe 

waves  which  converge  toward  0  This  suggests  that  on  the  surface  these  waves 

can  be  ropi'esented  by  Hankel  functions  which  r€pre''''n(.  re, n  verging  and  diverging 

waves  In  a  cylindrical  geometry.  This  behavior  is  also  suggested  by  the  represen- 

(1  2) 

tation  of  the  functions  E  '  {"U,  0)  In  eqn.  (4-14). 

m 

We  consider  the  tern  of  eqns.  (4-26) 


Ika(2ffif0-y2) 


lka(27ri+  ^  -0) 
ie 

^  sin  0 


(4-29; 
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This  la  a  valid  representation  awaj'  from  the  caustic,  for  ka(7r-0)'?';' 1,  and  la 
suggestive  of  the  asymptotic  form  of  the  Hankel  functions.  We  will  find  such  a 
Ilankel  function  representation  and  then  continue  eqn.  (4-20)  to  the  caustic  region. 
From  the  asymptotic  behavior  of  the  Hankel  functions 


H 


(1,  2) 

1 


+  l(z  -  ^^4  ) 


(4-30) 


wo  have  that 


•*■<?  -’^^2)  ^  -3t1/4  fj^  lka(2iri+  72  ) 

- -  =Ie  ,  2  ® 


and 


Jsln  0 


lka(2jrJt+  ^ 


It -6 


(4-31) 


[kv/v-oj] 


sin  B 


jsln  0 


If  now  v'e  substitute  eqtis.  (4-31)  and  (4-32)  In  (4-29)  we  get 


(4-32) 


-ljr/4 


Hour 
\  2 


ika(2jri+  ^2) 


IZ3 
'I  sin  9 


[ka(ff-O)]  g(?^KHj^^'[ka(jr-9)j  g(?'|)j 


(4-33) 
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Nov/  If  we  let  0  approach  t  we  have 


(4-34) 


and  (4-33)  becomcH 


u  ...  --  I  Ika(2»'-^  +  ^l2  ) 

Q  '  jkajJTs  c  2J\(kaOr-03  )g(?^)  (4-35) 

which  is  finite. 

In  the  above  treatment  we  have  performed  the  continuation  into  the  shadow 
caustic  using  the  g(^)  as  an  example.  Of  course,  the  same  will  hold  using  f(C). 

Using  our  new  representations,  eqn.  (4-33),  and  the  analogous  one  for  the 
f(5)'s  eqns.,  (4-2G)  become 

>  , 

I  ^  -3ri/4, - ITT  J  iM25ri+’^/2)r 

!  H^=H^sln0e  e  )f(?^) 


e  o 


+H{’^’(kaC»r-0j)f(5’^)| 


(4.-38) 


cos  0  e 


'/2  )|H®W|>-0])g(5j) 


where  as  above, 


+  Hj^^kHO-83 )  g(e;j) 

. 


and  8>r/2. 


m  ■  (ka/2)^/^ 

5^  =m(2;r.f  B  ~\) 
sm(2jri  +•—  -  0), 
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A  similar  representation  cun  be  found  for  the  higher  order  terms  near  the 
caustic  0=0.  We  give  the  results  without  repeating  the  analysis. 


„cw  „  ,  M  ffl/4  I - I  0  ^ 

II.  Sin  0  c  JjrrnJ—; — ~ 

0  o  ^  ’  Isin  0  Q 


6  Ika(27r?+  '2) 

(-)  e 


•  (kaO)  f(?^  ) 


=  H  cos  0  e 

p  o 


-1^4  3/2  p,  rT“ 

■"  ^  JHITe  —  <-> ' 


(4-37) 


•  ^  (kaO)g(f^^j)  +  Hj^^'(kaO)g(§^) 


where  0ijr/2  andka8ind<  1. 
4.2  The  Scattered  Field 


Again  with  the  field 


_i  _  A  -ikz 
E  =  E  X  e 
“  0 


(4-38) 


Incident  on  a  perfectly  reflecting  sphere  of  radius  a,  the  scattered  field  to  order  1/r 
is  from  eqns.  (2-16),  (2-17)  and  (2-25) 
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Wo  now  restrict  our  detailed  treatment  to  the  first  of  eqns.  (4-40)  and  con¬ 
sider  tl»e  sum 


e"  = 


Z: 


2V 


-1 


W2  ^  "■  y^(cos  a) 


5.,  1,  (ka) 


V-V, 


sin  ot 


(4-41) 


Eqn.  (4-41)has  the  contour  Integral  representation 

-1 

r  >//  .  (OA)  ;:)T’ 

2vdV 


■/ 


1+  e 


21Vjr 


^  l/(ka)  3P  ,  (cosa)  <p  ,  AV&)  P,,  ,, 

l/-y2  V-^l2  tJ- 12  ^-12 

(1) 


(cos  a) 


(ka) 


vr-  '^-1/2““’ 


9"  (ka)  Bin® 

(4-42) 

where  the  contour  C  encircles  the  positive  real  n'’ic  in  a  clockwise  direction  as  in 


figure  4-2.  As  before  we  note  that  the  contribution 


of  the  poles  at  v'=  1/2  vanishes. 


Since  2i!j^  -•  +  5^^  cqn.  (4-42)  can  be  written  as  tl»e  sum  of  two  terms. 


The  first  we  consider  if 

-1 


i^dV 
2Tvf 


1+e 


r  ^P P.  \/  (cob  a) 


V-V2 


3» 


sin  a 


00 


(nf^) 


,  fSP  \cOS  o)  P  \c08  O') 
IJ  n  .  n 


50 


+• — ~  “1 

sin  a 


(4-43) 


QA 
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00 

>^(n  +  ^  0-)=  I  cot  I  , 


hence. 


n=o 


.  d ,1  .ot .  ^  j. 

h'  da  ^  2^°^  2^  2  sin  a 


i.  “ 

1  i  =, 


(4-44) 


(4-45) 


From  this  resultjeqn.  (4-45)j  .ve  confine  our  attention  to 


e^/ 


VdV 


1+e' 


2lV7r 


L  ■>'-72 


da 


r<2)  X 

^-1/2 


p-l 

^-V2 


(cos  a) 


sin  O' 


(4-46) 


Since  the  integrand  In  eqn.  (4-46)  is  an  odd  'notion  of  V  we  can  reflect  the 
lower  part  of  c  in  the  origin  and  get 

(2) 


e 


oofie 

s  / 


-urfit 


2Va-i 


(2)' 

(ka) 

L  V'  V2 


9P  \  (cos  a) 

^-V2 

da 


^'-V2 

?  ,  (ka) 
1^-72 


P  !,  (cos  a) 

p-n 


sin  a 


(4-47) 


where  e  >  0  is  a  small  parameter.  In  the  sequel  drop  the  c  in  the  limits  of  inte¬ 
gration  with  the  understanding  that  the  contour  is  to  run  Just  above  the  real  axis. 

As  in  the  previous  discussion  of  the  fields  on  the  surface  of  the  sphere  w'c 
would  like  to  close  the  contour  in  eqn.  (4-47)  by  a  large  semi-circle  in  the  upper  hal^ 
nl'hie*  To  this  end  we  decompose  the  Legendre  function  as  in  eqni  (4-lu), 


-  2P“\  (cos  or)  =  a)  + 

^-1/2  1  1  '  ' 


(4-48) 
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and  requii'e  o^O.  Since  E  ^  {v,a)^  ...  e  is  the  dominant  term  in  the 
I-egendre  fiinction  in  the  upper-half  V-plane, 


2i  t/ir 


-2it/ IT 


1  rl  \ 

=  1  +  0(c  ), 


(4-49) 


and  since  in  parts  of  this  half  plane 


VV2  ^ 

-J3.-  =  0(1) 

5  ,  (ka) 

V- !  2 


(4-50) 


we  have  that  the  terms  of  the  integrand  of  eqn.  (4-4?)  are  exponentially  small  for 


1  except  for 


r  5®;.  (. 

1  '*'“72 


(to)  ' 


[  V-/2 


(11/  <>“> 

*'-72 


(2) 

e;  q) 
sin  a 


(4-51) 


The  remainder  is 


00 

■h 


_ _ 


1+e 


2t^jri 


(1) 

?  ,  (ka) 


sin  a  1  +  e^ 


E^j\v,a) 

sin  a 


in)  -1 1 

^E‘"(V,Qr)  j 

Bo  ..  -2V5ri 

J 


(4-52) 


1-t-e  j 
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Eqn.  (4-52)  can  be  determined  by  closing  the  contour  and  evaluating  by  means 
of  the  residues  at  the  zero  of  (l<a)  and  ^  (iai)  .  if  we  designate  the  zeros  of 

(1)  "-V2  ’'-V2 

,  (ka)  by  V 

Jl) 


(i)' 


and  those  of  ?  (ka)  by  1/ 
V  -1m  1 


72 


m 


,  (ki.)=0, 


then 


.<2). 


sin  a 


1 

1+e  ” 


Ef’(Va> 


Sin  a 


^  (1) 

.«) 

i  m 


1  +  e 


+  Tr\ 


(2)' 

5;  1/  (ka) 


m 


■"fj.  5",  to)] 

2/.I/2 


30 


1+e 


,0) 

1  m 

3o 


i+e 


2y^Ti 


(4-53) 
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Since  ,  and  its  derivative  vanish  at  the  poles  we  can  replace  in  eqn.  (4-52) 
V-HO  y-V2 

by  .  Also,  since  the  zeros  of  and  start  atll^l/'-'ka  for  ka  snfflc- 

U-ri  *'-1/2  s^-V2 

iently  large  we  can  use  the  representations 


~  "  60  “  ^  [t^v'(t)  +  4t  v(t)j| 


(4-54) 


where 


t  =  —  ( V  -ka),  m  ^  (ka/ 2)^^  ^ 


land  v(t'  and  •'.Vj(t)  are  the  Airy  tunctions 


v(t)  =  fr  Aj(t), 

wj(t)=,nr[Bj(t)+iAj(t)]  . 


Since  Im  there  are  the  convergent  e^qiansions 


—  -T<^‘ 

2TiV 

^  i=0 


23ri1'.? 


1 

-23ri3^ 


^=0 


i  2jrlVi 
‘  e 


(4-55) 


(4-56) 
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From  equations  (4-17  and  (4-18) 


{-ii.  a) 


00 


2wi1^ 

1+®  £=0 


Ef  ^  Q) 


00 


(4-57) 


1+e 


-2i[iV 


=  i  ^ (t^.  2tC^+1]  -q) 


and 


£<J^(1/,  a)  ^  - 


i=0 


r<^-b 


P  (Vt  1  )  ij  r  sin  a 


i(i;  a  -  ) 


(4-58) 


So  on  substituting  eqtis.  (4-54),  (4-57),  and  (4-58)  in  eqn.  (4-52) 


on 


,  ,2  i 

I  -{&'-)  1/2  ® 


Vi/4 


m 


(4-59) 


where  we  write,  noting  that  this  is  a  different  notation  than  the  previous  for  ^ , 

^  =  (2xJl  +  T  -  0)  m  = 


^  -  (2rJl  +  X  +  0)  m  =  wip^ 
/•“  Jft 

^  j?" 

-00 


(4-60) 


wj(t) 


dt 


Here  we  use  the  caret  notation  not  to  designate  a  unit  vector  but  to  be  consistent 
with  the  notation  of  Logan  ^95^  . 
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We  return  to  eqn.  <4-5l)  and  make  a  saddle  point  evaluation.  The  staC.onary 
point  occurs  at 


ka  cos  0-/2  or  ka  sin  0/2 

So  by  the  standard  "lathods 

P  j  -2ika  cos  0/2 
l"  =  -kae 

Finally  we  have  the  egression  for  the  fields 

e“^fl  -21kaco8®/2  (ka)^ 

E  =  E^icos0-i^^jkae  +-1/2® 

I  ID 


i=0 


E  isinp  ® 

1^0 


ikr 


kr 


[f 


,  -2ika  cos  %  3ri/4 

— ~^f7.  ® 


00 


)  e  I 

je=o 


[where  is  evaluated  analogously  to  £|  and  Ub  functi(H>  (><€)  is  given  by 


00 


P(5)  =  l/fF 


/ 


-00 


i?t  v(t) 


Wj(t) 


(U 


(4.61) 


(4-62) 


(4-63) 
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Near  the  forward  direction,  ,  it  is  necessary  to  modify  our  treahnetrt 

(1  2) 

because  of  the  singular  behavior  of  the  functions  E  j  ’  o).  We  will  start  with 

the  integral  representation  of  eqn.  (4-42)  but  retain  the  Legendre  hmcticna. 


^  t  2vAnl  VV2^^' 

5.  %a)  sino 

'c  L  1/-V2  1^-72 

where  c  is  a  clockwise  ccmtour  about  the  positive  real  axis.  On  reflection  of  rhe 
odd  part  of  the  lower  contour  this  becomes 

ka  [5^^^  (ka)  3P~\  (coso)  (cos  or)  1 


f*"'  (W 
.  I'-Va 


r  1 ,  V 

>^-72 

90 


«^-V2 
5%  (ka) 

H-H2 


t^-V2 

sin  a 


I  OQ 

i  i 


2vd1>\ 


!  ***  ^  /tv 


L'v-y2^‘“' 


0) 

5  ,,  (ka) 


P  \  (cos  or) 

i^-V2 


dP  \  (cos  a) 

V-V2 


p”l,  (cos  or) 

y-y2 


-a>H€ 


1  to 

P  ^  (cos 

V-V2 

■^-72 

da 

5  ,,  (ka) 

v-Vz 

sin  or 

(continuedl 
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(cos  a)  p”l  (cos  or) 

^-y2 


2jrW[  Jor 


vdv  3  (cos  or) 
l+g“2ril/  _  ^of  sin  O' 


(4-65) 


The  first  two  terms  in  eqn.  (4-65)  are  evaluated  using  file  Bessel  fiicctioa 
representatimi  of  the  Legendre  fUnctimi 


P  ^  (cos  O’)  =  5  J  <  vq) 

^‘V2 


(4-66) 


for  \V  {  large  and  a  small.  The  result  is  for  the  first  two  terms 


ka  5’^!,(ka)  p'^  (cos or)  *a>  JV',  i/  (to 

V-»  t“' »)  4.  (\«-) 

V-k  ^-/2 


y.®  2V',  If  (ka)  p“^  (cos  or) 

/  ^-^2  V-V2  _ 

/  (1)  " 


J,(Qfcl) 


=  -  2  ji*  ■  m  jika  p(0)  +  m  p'^'  (0)  +  ...  I  (4-67) 


Ivdiere  we  write 


P(?)  = 


2i  ^ 


0  o 

J-  /  i«  „  1  f 

'  /  "iW  f'  J 


*Pt  v(t) 

e"  — --  dt  (4-68) 

w,(t) 


(1)  3p(?) 
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aiidp®(E)= 

’v 

S'-Vs 


Similarly, 

- if - »dy+  /  — —  — if - 2v&tf 

J  . 


=  2{x~ j|  (ka  a)  m  q(C)  +  m  q^^^(O)  +  . .  .J  (4-69) 


viiere  the  q's  3U*e  the  functions 

0 


^  F  J  Wj(t)  iTTJ 


JD.eX- 

q  (?)  -  gg 


(4-70) 


j  The  third  tc;rm  ut  jqn.  (4-65)  is  evaluated  by  using  the  small  ang^e  e^qiansionJ 


P  y  (cos  a)  =  1  -{V^-  2  )su?  ^  + 
V-n  4  2 


,  .2  1  ,2  9  ,  .  4/r 

(V  -^)(V  --)sin  - 


+....  (4-71) 


IWe  find 


r  y-V2  V-V2 


-  ^  (ka)^-(^  (ka)^-  |(la)^sin‘‘ 


•*■(24 


(4-72) 
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The  fourth  term  is  eval:iateH  using  flie  Bessel  hmction  representation  of  the 

i 

Legendre  function  and  expanding  tiie  denominator  (l+e  “  ).  Again  we  remailc  tiiat 

(2) 

integrals  of  this  form  can  be  evaluated  by  residues  and  hence  we  can  r^lace  <  ^ 

h 

by  2^  SO  as  to  keep  the  form  standard.  Perf<»ming  Oiese  opci'ations  and  ke^iz^ 

V-V2 

the  hi^er  order  terms 

Iks)  9P  ^ ,  (cos  or)  (ka)  P  *  (cos  or)  1 


g-2kajlFin y  j|(ka«)  4[m(l+l)2r] 

i»0 


where,  as  before 


a 


m  v(t)  ^ 
B  ,lx  at 


w|<« 


The  last  two  terms  d  eqn.  (4-63)  are  equal  and  can  be  evaluated  usiiig  the 
ie»parsi<»  in  e^  (4-71).  We  get 
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5P  I  (cos  a)  p~l  {cos  a) 

■if- 2  y-Va 

-  ^^41  _ — _  .  _ -  _ 


/dv  ,,  icosflr; 

1  y-  72 

- - -  -vdV - + 

i^-2Tyi  JO 

u 

J  C 


7 


(cos  or)  P  ^1.  icosor) 

2 _  ^  )  1^-72 


1  sln^  -SL+  JgL?6?- 

12  170  2  4,193,530  2  ’ 


rhe  fields  are  given  by 


ikr 

e“=  E  ico8  0-^ 
d  o  kr 


|^(ka)2  -  f  siii^|+ ... 


J,(oto) 


+  2^jr  -  m  |lkap(0)  +  in  i»  *«-•>  +  ...  j 

oka  1-  .if 

+  2^jJ(Qka)m  [ika <i(0)  +  m q^No)  +  . . . J 


2  ka  fr~  m 


(4-75) 
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£ 

o 


J.{rfca)  - 

+  2(^  m[ikaq(0)  +  mq'^{0)  +  ...j 


+  2Jt  j|(aka)m  [ika  p(0) -i- m  p^No)  +  . .  .J 


00 


-  2ka  fw  m  ^2T(f +l)k»  ^  |^i(f +l)j 

/=0  ^  " 


(4-76) 


4,3  The  FoniMilas<rfSectlTi!s4.1aiid4.2. 
For  a  plaiiC  elo'^tnanagnetic  wave 


E*  =E  xe"**®  (4-77) 

“  o 

incident  on  a  perfectly  condactisg  q^re  of  radius  a  we  present  tiie  formitlas  of  the 
previous  sobsectiaas  along  with  a  brief  commet  t  on  tbe  physical  interp.-etattoo  and 
the  methods  calculattoo, 

4.3.1  The  field  on  the  sur&ce 

I  In  the  lit  region  including  the  caustic  we  have  found  d>e  fielii  t?  oon^t  of  cwo 
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terms.  The  first  is  that  due  to  direct  illumination.  This  is  diaracteiized  by  the 
fact  that  (1)  it  remices  to  gecnnetric  optics  in  the  limit  ka — and  (2)  it  carries 
the  i^ase  of  the  incideiit  field.  Under  die  smnearbat  arbitrary  eonditioa  0 ^  ^60° 


we  have  the  direct  field 


H^=-2H  c-3e©sin0e~*^®®*® 


{is! 

1+ - 

2kac 


ain^  e  5  sin^e  -sin^d 


cos  $  2(ia)  cos  9 


H?*=-2H  coe0e 

f  o 


-ikaeosd 


i  sin^9 
2kaeos^9 


9  sin^9-  8in^9 


20ar  cos®© 


(4-78)i 


I  Hie  second  contribotioo  in  this  regioo  arises '.'em  wa^s  that  have  cr^ 

I 

around  the  bade  of  the  qibere  and  hence.  Is  cfaaracterixed  by  having;  the  phases 
3t 

Ka(2si-i‘'^  >f(»r9>0.  The  fmm  we  give  £<»■  fids  contribittioo  depends  ^>on  the 
'-sine  of  ka  sin  9  althoogh  we  iznpwe  the  restrictiop  9^s/2. 

For  6^-k/2  Cos  a,  wave  coatriiaitiOB  is  for  ka  sin  9  >’  1, 


o» 

i-0  I 


|"=-H^cos0 


MS  hilt 


(4-7S/  I 
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and  for  ka  sin  9  <  1 


cw 

H^  =  e  H  sin 
B  o 


ika(2ai+  ^^2) 


•  |Hj^^'(lae>  Oafi)  UsJ  )  j 

«  -i'/^  3/2  T,  J 

l^=e  H^cos^m  jrj— ^  {-/e 


(4-80) 


and  fffi*  0  =  0 


*  (kaO)  g(f^+jH  (kaO)  g(5j 


,  ir/4, _  r’  j  ikM(2ri+^m 

(-^'c 


(4-81) 


H^=2H 

f 


-i»/4  ,/2  I  Ua{2ai+^/2)  , 

[^cos|»e  •  {-  )*o  fsaJi  g(€,^,) 

^=0 


In  eqns.  (4-80)  and  (4-81)  tiie  bmctioa  f(S)  and  g(C)  appear,  Tt^u:  have  been 


coo^iuved  and  tabolated  by  N.  A  Logan  ^ISosj 


In  the  transition  region  30°i  0  i90®  there  are  again  two  cootrarailcns.  t  be 


first  is  the  c<»tinuatioo  of  the  shadow  currents,  file  seconi  is  fiie  crewing  waTe 
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contributioii.  These  are 


-Dca  co«  0  „ 
e  i(-ni  cos  6) 


0  ^  “'■* 
^  1  -ll(»coe« 

H4=-H  cos^r===  e  f(-mcoee) 

0  o 


(4-82) 


CO 


Z.V[ 


i«o 


Here  we  have  made  the  aubstitatton 


"^+1  ^  1 


5^  =  m(9  - 1  =  m  sin  (8  - 1 ) 


=  -  m  cos  0 


In  the  first  term  so  that  fids  term  has  the  pk'se  ct  the  incideBt  wave.  The  sfibsequei*. 
terms  are  just  fi:e  creeping  ^tc  cc:itributisas  realized  fimn  their  f^ases. 

In  the  ^ladow  region  away  fixmi  fite  caustic,  8 ^  r/2  ,  ka  stn  8  >  1. 
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These  are  purely  creepii^  save  tenns  agalc  cfatracterized  bv  <Jieir  pluses. 
Finally  near  tite  slsdov  cau^c  «e  have  for  ka  sin  0  ^  i,  0  >t/2. 


<4-83> 


iica(2»j(+^J 


[ka<»-0)]  )j 


(4-84) 


„ _ jwl/A  3/2  _l  a-0  ^..1  ' 

fi^=co.0H^e  m  r\^ ^  <-^  « 

/=0 


I  iia{2»/+%) 


a»0  =  » 


H  *2H  sin 
9  O 


H^-2H  c<» 

f  o 


<4-8S) 


^.5/a  ,  ika(2»i+%)  _ 

{-4^  Jj[la(»-8)]g(5^) 
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4.3.2  Hie  Far  Field 

In  tile  &r  field  tiiere  are  essoitiiily  just  two  legioes.  Hie  first  is  tiie  large 
r^<m  for  which  The  second  is  tiie  f<nward  scattered  region  far- 

G'-i  and  ka  sin  1. 

For  0^5  the  dotuinant  contribiUlon  will  arise  frem  tiie  terms 


-2ika  cos 


')f^ 


_  .  2e 

i  sin  — 


Wcoe«f 


(4-36) 


''  8\  ^e^/ka  -2Ikacoe%\/ 


2kacos''| 


2  0 
V 

iHdcfa  in  tLe  limit  ka— reduce  to  the  geometric  0|)tics  fields.  Tl»  adifitioial 
terms  c<»rres|>and  tc  tiie  creeping  wave  comribution  and  are 


ikr 


(4-6i) 
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where  m  =  (ka/2) 


=(2B-i+  r-0)  m  =nn^ 
=  (2t-/+  r+d)  m  = 


-00 

■^f 


iSt  v(t) 
e  ZITT  dt 


i5t  ^ 


arc  functions  which  have  been  extensively  tabuiaied  by  Logan  The  creeping 

wave  contributions  are  a  very  small  correction  except  for  cases  in  which  there  is  a 
phase  correlation  between  the  primed  and  u]q>rimed  terms  in  the  brackets.  This 
occurs  for  jr  or  mod  2jr.  In  facv  5or  0-0,  baclwcattering,  this 

correlation  approach  can  be  used  to  predict  the  relative  maxima  and  minima  in  the 
backscatiering  cross  section  as  a  function  of  ka. 

In  the  case  of  forward  scattering  there  aiv.  .Jjc-  ~'^o  distinct  contributions.  The 
first  can  be  recognized  if  we  recall  the  cross  section  theorem  which  states  that 

or-  =  -^  Im  ffx) 

where  is  the  total  cross  section  and  f(0)  is  the  complex  field  aii4>iitudG.  Since  in 
jthe  limit  of  geometrical  optics 


cr^  =4xa 
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we  have  in  this  limit  that 

Im  f(a-)  ={ka)^. 

Recalling  the  form  of  the  fonrard  scattered  field  we  have  the  optics  contribution 


=  iE^cos0^j^  j  (ka)^ -Q(ka)^-|(ka)^]  sin^or/a)  +  . . .1  (4-88) 

\  "  >t)p  J 

is  the  same  eiiqiression  with  cos  0  replaced  by  sin  0 


an( 


^E^^  =iE^sin0^j^  |(ka)^-  Q  (ka)^- 1 (ka)“]  sin^(tt/2)+ ...| 


(4-89) 


The  remaining  terms  in  the  forward  scattered  field  are  those  that  arise 
directly  from  the  shadow  boundary 

r.T  trAraS 

p(0)  +  mp'*\o)  +  . . .  J 


ikr 


/  »\  fjJaJta)  _  ... 

(  E®  J=  i  E^cos  0^-  (2 /Fm) -  — '  ’• 


aka 


j|  (a  ka)  |ika  q(0)  +  (0)J^ 


(4-90) 


Jj(a  ka) 


I  E  j  )  =  i  E  sin  0  (2  jjw  m)\ - - —  P-m.  qlf;)  -r  mq'^\o)  +  . .  .1 

N.  ^/'SB  aka  J 

+  jj  (a  ka)  juia p(0)  +  mp^^\o)j  | 
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V 

THE  PHYSICAL  OPTICS  APPROACH 


Perhaps  the  best  known  technique  for  obtaining  approximate  solutions  of 
high  frequency  diffraction  problens  is  the  method  of  i^ysical  optics.  The  key 
feature  is  an  assumption  about  the  current  distribution  on  the  surface  of  the  scat¬ 
tering  object,  and  in  this  section  the  mt:thod  is  applied  to  the  case  of  a  perfecay 
conducting  sphere  cf  I'adius  a.  The  degree  of  approximation  involved  is  examined 
by  comparison  with  the  exact  AKe  series,  and  a  numerical  example  is  treamd  which 
lends  support  to  the  use  of  {diysical  optics,  particular  where  the  main  purpose  is  to 
obtain  general  estimates  of  the  scattering  behavior. 


5.1  Physical  Optics  for  the  Stohere. 

The  scattered  magnet'.c  field  Is  given  In  terms  of  die  current  J  induced  in 
the  surface  of  the  sphere  by  the  equatitm 


I 

i 


(5-1) 


where  R  is  a  unit  vector  from  the  receiver  to  a  v?-  '  ’.bJe  point  (a,  G',  0' )  <m  die 
^here,  the  distance  between  these  points  being  denoted  by  R.  If^  were  accurately 
known,  die  above  equatimi  would  provide  an  exact  expresslmi  for  the  scattered  field 
and  die  basis  of  the  physical  optics  approach  is  an  si^roxinmiiim  to  the  true  value 
ofa 
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According  to  ray  theory  thei-e  exists  a  sharply-defined  shadow  regi!.  n  bshind 
the  sphere  in  which  the  total  field  is  zero,  and  since 

J  =n  AH  > 

where  ^  is  a  unit  vector  normal,  the  current  distribution  over  the  shadow  area  must 
be  identically  zero.  For  the  illuminated  portion  of  the  inhere  J  is  obtained  ou  ihe 
assumption  diat  the  field  is  reflected  at  every  point  as  though  an  infinite  plane  wave 
were  inclrlont  on  as  infinite  tangent  plane,  and  this  givea 

J  =  2n  -N  (5-2) 

that  is,  twice  the  Uiagential  component  of  the  incident  magnetic  field.  The  current 
distribution  is  now  completely  specified  by  dioosing  the  incident  field  as  that  given 
I  in  eqn.  (2-1),  and  hence,  at  a  point  (a,  8',  0’)  on  the  Illuminated  side  of  the  qphere, 

I  J  =  (cos  8’  X  -  sin  0’  cos  0'  z)  2  H^e  ^  ®  .  (5-3) 

Moreover,  the  fact  that  r  is  large  compared  with  the  radius  of  the  ^here  means 

ft 

that  R  i."  .  ffectively  directed  toward  the  origin  and 

R  r  -  a  cos  9  cos  8'  -  a  sin  8  sin  8*  cos  (0  -  0’) 
from  \diidi  it  follows  fiiat 


1 
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7/2  2r 


f  f  .  ^ 

=  j  I  sin^e’  cos  0'  d0»  dfl* 


0  0 
7/2  2» 


..// 

0  0 


e"^  sin  e’  cos  0*  d0’  d0’. 


The  component  is  of  the  order  l/r  and  dierefore  negligible  by  comparison. 

The  above  integrals  can  <xily  be  evaluated  exactly  in  certain  qpedal  cases 
and  for  arbitrary  values  of  6  and  0  it  is  necessary  to  rely  up<m  approximate  tedi- 
niques  based  upon  die  (assumed)  larj^  value:  of  ^  These  tediniques  can  be 
illustrated  by  reference  to  I^.  Here  the  0*integral  is 


27 

f  ^-ika  sin®  sine'  «.s(0-0')^ 

2t-0 

/-ika  sin®  sin®'  < 


^  ika  sin®  sin®  cos  0  sin0*)  d0’. 


and  since  die  term  involving  sin  0'  contributes  nothing,  the  integralion  being  over  a 

conqilete  period,  we  are  left  widi 
27  "0 

4  t  -ika  sin  ®  sin®'  cos  0'  «  ,  a  »/._>„  1  «t\ 

cos  V  f  ^  ^  cos  0  <*0  =  ~27i  cos  0  Jj(ka  sin  0  sin  ®') 
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and  hence 


f 

■  -2r.i  cos  j}  / 


-ika(eos  ^l)cos  S’  n  -  a  » 
e  Jj(lca  sm  S  sin  S')  sin  S'  dS'. 


If  it  is  now  assumed  that  ka  sin  S  sin  S'  is  large  compared  with  uni^'"  the  Bessel 


functicm  can  be  r^lacsd  oy  its  asymptotic  expansion  to  give 
_ _ 


L  -2i -  cos  0 

^  jlkasinS 


-ikafcOS  SM)C0S  S' _ .  a  ^  o,  3x  .  , 

e  c<w(ka  sinS  sinS'-  —  )sin  S'dS' 

4 


and  bv  writing  ooh  {ka  sin  S  sinS'  -  ejmonentie!  form  it  can  be  verified  fi>at 

4 

die  raiy  saddle  point  for  whidi  0  jf  Six  and  Oss  ^  is  S'  =  S/2.  This  is  equivale: 


to  substituting 


—  exp  i  [•^  -  ka  sin  S  sin  S'j 


4  *3^  the  cosine  fiictor  and  gives  ^  q 

I  ix/4  r  2x  a  ^  ^  -2iks  cos  7  cos  S*  .  3/2,  S  . 


1  s’  ^  ® 
{tocos-  Ig-J) 


^  S_  .  -2ikacos| 
B  ^  tan ~ cos  we  * 

ka  2 


e“'dt 


coa^  2 


The  fidlure  of  this  cooditioti  at  the  lower  limit  of  IntegratioD  clearly  indicates 
that  even  U  fim  subsequent  evaluation  of  1^  and  I2  were  performed  exactly  die 
results  would  at  best  be  approximate. 
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Providing  B  r  (cases  alreatfy  excluded  by  the  requirement  that  ka  sin  9  sin  0* 

be  large),  the  limits  of  integratton  can  be  replaced  by  +  a>  and  since 

2 

f 

-00 


we  finally  obtain 


A 

-2iJca  cos  -j 


^  9  ^ 


An  analogous  treatment  applied  to  fiie  integral  leads  to  the  result 


T  ly 
*2"  ka  ® 


9 

•zusa.  cos 


and  if  these  expressions  for  Ij  and  are  inserted  into  eqns.  (5-4)  and  (5-5),  the 
scattered  field  takes  the  form 


a 

O 


cos  0  e 


-2ika  cos 


9 

2 


(5-6) 


Ikr 


_8  a  _  e ,  ^  -2ika  cos  o 

=  w  E  - sin  ^  e  ^  . 

p  ^  o  r 

The  corre^KMiding  scattering  crctts  actions  are 


{5-V) 


and 


2s 

CM  P 

8in^^ . 


(5-8) 

(5-9) 
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The  above  results  have  been  obtained  by  an  approximate  evaluation  <d 
and  1  ,  the  basis  of  wMdi  is  the  assumption  that  ka  sin  6  is  large  ccm^pared  witti 
unity.  The  smaller  the  value  of  sin  6,  the  larger  must  ka  be  in  order  to  hilfiU 
this  requirement  (sr.o  footnote  p.  1 19  )  and  in  the  limiting  c-.  es  for  which  sin  d=C 
the  method  is  not  Icmger  valid,  ft  is  fortunate  that  tiiese  cases  .'.re  the  very  ones 
for  which  Ij  and  can  be  treated  exactly. 

d  =  r  correr^Kmds  to  forward  scatter  and  since  ^  is  dien  zero 


I,=0 

A 

and 


giving 

and 

The  ptdarization  Is 
cross  aecUtms  are 

and 


_.s  ika^  „  e^  ^ 


2  ikr 

"oV  ■-“I'- 


U:  r  ::fcrc  identical  to  that  of  the  incident  £eld  and  the  scattering 


1.2  4  2^ 

o  =  T  k  a  cos  p 
9 

1.2  4  .  2^ 

O0  =  v’ka  sin^  , 


implying 


a  =  ag+a^  = 


4aA 
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where  A  =  7ra“  is  the  geometric  cross  sectional  area  of  the  sphere. 

The  other  limiting  case  is  0  =  0  and  correiponds  to  back  scatter.  Since 


^  =  2  cos  O', 


lj=0 


.  /  -2ika  cos  9'  ,  n  f  -2ika,  ^  -  -2il;a  . 

1  =27:  f  e  sin  0'  cos  0'dO'  =  -  - - <  e  + - (e  -i) 

2  /  ika  21ka 


if  ka  is  large.  Hence 


)r  -2ika 
—2—  e 

ika 


ikr 

a  „  e  -21ka 

_  cos 


_s  a  _  e  .  -2lKa 

E  j  =  -  E  -  sin  0  e 

0  2  o  r  ^ 

and  these  are  in  complete  agreement  with  equatio.is  (5-6)  and  (5-7)  notwithstanding 
the  fact  that  the  approximate  method  of  evaluatlni;  ana  breaks  down  when  04. 

In  view  of  this  continuity  as  0  approaches  0  it  is  reasonable  tc  put  forward  the 
scattering  cross  sec  tions  given  by  equations  (5-8)  and  (5-9)  as  valid  for  all  0  not 
near  to  tt  and  the  implications  of  such  a  statement  wlil  now  be  considered. 

2  A  Comparison  of  Formulae 


According  to  physical  optics  the  scattered  electric  field  at  any  point  consists 
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ol'  two  in-phase  rectangular  components  in  the  aperture  plane  of  a  receiver  directed 
towards  the  center  of  the  sphere.  This  is  certainly  in  partial  agreement  with  the 
conclusions  of  the  exact  analysis. 

As  regards  the  nature  of  the  scattered  field  and  its  dependence  on  jj,  the 
predictions  of  the  physical  optics  method  are  correct,  as  may  be  seen  from  a 
comparison  of  equations  (b-6)  and  (5-7)  with  (2-37).  But  the  approximate  treat¬ 
ment  has  produced  components  which  are  in  phase  for  all  6  and  moreover,  has 
destroyed  most  of  the  dependence  on  6.  Indeed  if  we  exclude  for  the  moment  the 
case  of  6  neai-  to  n  use  of  tne  physical  cities  methr.i  equivalent  to  replacing 


and 


V’  —T  ..-‘(cos  e)] 

(-1)""“  aJ-S - 

1  Z_J  n  L  sin  e  J 

(cos  0)1 

B  - > 

n  sin  0  J 


(5-10) 


n=i 


00 


f  9P*  (cos  0) 


(5-11) 


ka  -?-ka  cos  w 

by  X  -7  e  ^  respectively,  where  A  and  ^  are  defined  in  eqn.  (2-25). 

n  n 


In  particular,  for  backscatter  it  replaces 

CD 


.  „  ka  -2to 
by  --g-e 


n=l 


and  for  the  exceptional  case  of  forward  scatter 

n=l 


\jkj 


i‘ 
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The  degree  of  approximation  which  these  Imply  has  been  examined  for  a  selected 
value  of  ka  and  the  results  are  given  In  section  5.  3. 

Such  a  comparison  Is  not,  however,  a  fair  test  of  the  physical  optics  method 
In  that  additional  appi  oxlmatlons  were  made  In  order  to  evaluate  the  integrals  and 
The  basic  assumptions  as  to  the  form  of  the  currents  can  only  be  tested  b>  a 
study  of  the  currents  themselves  and  this  will  now  be  done. 

At  any  point  (a,  6,  0  )  on  the  surface  of  the  sphere  the  physical  optics 
approximations  to  the  current  can  be  obtained  from  cqn.  (5-3)  as 


.1  =0 

(5-12) 

r 

T  „  „  -Ika  cos  0 

J.  =  2  H  cos  0  e 
d  0 

(5-13) 

and 

r  n  X,  t  A  a  “ika  COS  0 

Ju--2H  sin  0  cos  $  e 
p  o 

(5-14) 

for  0^0  ^  jr/2  (illuminated  portion  of  sphere),  with  J  =  J-  =  =  0  otherwise. 

rap 

In  contrast  the  exact  current  dJstribnlion  is 


J=n>^(H^  +  H®)  (5-15) 

i  s 

where  H  and  H  are  given  by  eqns.  (2-181  (2-22)  aud  (2-25).  When  substituted  Into 
eqn.  (5-15)  these  e^qjresslons,  together  with  the  Wronskian  relations 

[ka  J  (ka)  j  l/ka 

j  (ka)  -  h  (ka)  /- - ^ —  =  - 

n  n  ^kah^(ka)j'  [ka  h^(ka)] 
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j  (ka) 


lead  to  the  formulae 
00 


3P^  (cos  0) 

2nrt  I  ^  ,  _3f _ 

n(iH-l)  ain  0  ^kah  (ka)j' 


(5-16) 


-j^  sin  ft  2^(- 


dP  (cos  0) 

r  1  p 

P  (cos  0)  QQ 

ir^  ^  t— ^-1 - ‘ _ 

IL  11  J  n 


(5-17) 


A  comparison  of  equations  (5-13)  and  (5-14)  with  (5-16)  and  (5-17)  now  shows 

that  the  physical  optics  approximation  to  the  current  replaces  ^ 

oo  9Pp(cos  0)  ^ 

_  1  y"...n2nH  f  ,,  _ 1 


a=i 


n(tt4-l) 


+  i  r- 


ka  h  (ka)  sin  0 
n 


T  •;  - 

2  ka 


.iH-l  2n'fl 
n(nfl) 


[ka  h_^{ka)]’  J 


P  (cos  0) 
n 


^^(cos  0) 


fka  h  (ka)l '  sin  0  ka  h  (ka) 
u  n  J  n 


(5-18) 


(5-i:;; 
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I„  -ika  cos  j  „  „  -Iks.  cos  0 

2e  and  -2  cos  6  e 

respectively  for  Q^6^ir/2  and  by  zero  for  other  values  of  6. 

It  must  be  emohasized,  however,  that  the  usefulness  rC  be  physical  optics 
approach  in  scattering  problems  does  not  depend  entirely  upon  the  accuracy  of  the 
current  distribution  which  it  predicts.  The  parameters  of  most  practical  impor¬ 
tance  are  the  far  field  amplitudes  and  the  fact  that  these  can  be  e;q>ressed  as 
stationary  forms  involving  the  currents  (as  in  the  variational  fomiJation)  suggests 
that  sli^t  errors  in  these  currents  do  not  necessarily  reveal  themselves  as  errors 
in  the  far  field  amplitudes.  Ideally  it  would  be  desl’-ible  to  carry  out  a  direct  com¬ 
parison  of  eqns.  (5-4)  and  (5-7)  with  (2-37)  with  no  approximations  made  to  Xj  and 

I 

but  the  labor  involved  in  a  numerical  integration  of  these  integrals  prohibits 

such  an  underiaking.  The  currant  distribution  is  the  only  altemstive  basis  of 
i 

}  comparison  not  involv  ng  approximations  additional  to  those  of  the  physical  optics 
method  itself. 

5. 3  A  Particular  Case 

A  significant  test  of  the  predictions  of  physical  optics  can  be  achieved  by 
confining  attention  to  a  single,  judiciously-chosen  value  of  ka.  The  case  wh^ie 
ka  =10  is  convenient  for  computational  purposes  aiid,  in  addition,  leads  to  a  b:^V,: 
saa  tcering  cross  section  whose  exact  value  (see,  for  example.  Key,  Scewart,  Pinson 
and  Prince,  [lOS^  )  differs  from  »hat  of  physical  optics  by  (about)  the  local  mean  of 
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these  deviations  as  a  function  of  ka.  To  this  extent  it  is  representative.  Woreover, 
it  corresponds  to  a  sphere  of  sufficient  size  for  the  results  to  be  of  practical  use, 
but  small  enough  {a  r::  1. 6^  to  give  a  stringent  test  of  the  physical  optics  method. 

The  basic  assumption  of  physical  optics  is  an  approximation  to  the  current 
distribution  which  would  seem  to  be  justified  if  all  dimensions  oi  the  body  (ino^dlng 
the  radii  of  curvature)  are  large  compared  wilh  a  wavelength.  Neverthelesn,  the 
method  is  known  to  give  good  results  for  a  wide  variety  of  bodies  not  excluding 
those  having  point  singularities  or  sharply  curved  surfaces,  and  indeed,  a  sphere  of 
radius  falls  into  the  latter  category. 

A  comparison  of  the  postulated  currents  witu  their  exact  counterparts  for 
j  such  a  sphere  (see  figures  5-1  and  5-2)  reveals  a  remarkable  amount  of  agreement 

f 

I  over  the  entire  illuminated  surface,  the  only  real  discrepancy  being  near  to  the 
!  shadow  edge  in  that  current  which  is  assumed  to  be  discontinuous  thera  Over  the 
shadov/  area  the  curreut ,  are  not  zero,  contrary  to  assumption,  but  the  amplitudes 
are  appreciibly  less  than  for  tlie  other  hemisphere,  particularly  in  the  case  of  the 
'continuous'  current  J^.  At  0  =  x  the  airrents  are  identical,  with  their  amplitudes 
showing  a  mai'ked  increase  as  this  point  is  approached. 

The  failure  of  tlie  physical  optics  approximation  to  the  current  in  the  snadow 
IS  not  surprising  since  the  currents  here  have  to  'fit  in'  with  the  unnatural  foi-m 
I  f  sroed  upon  them  in  the  other  region.  Moreover,  the  discrepancies  are  unlikely  to 
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exact  analysis 


FIGURE  5-1;  AMPLITUDE  OF  CURRENT  J-  (ka  =  10) 


T  1  exact  analysis 


29 


FIGURE  5-2:  AMPUPTUDE  OF  CURRENT  J, 
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have  much  effect  on  near-backscatterlng  In  view  of  the  shielding  action  of  the 
illuminated  hemisphere. 

The  fact  that  the  overall  agreement  between  the  current  distributi«His  is 
greater  than  had  been  s'jpposed  does  suggest  that  the  discrepancies  >«4iich  have 
been  found  in  physical  optics  values  of  scattering  cross  sectio:'.c  are  not  necessarily 
attributable  to  errors  in  the  currents  themselves.  It  may  well  be  that  for  bodies 
having  no  surface  singularities,  and  with  a  receiver  in  the  illuminated  half-*^acc. 
the  major  Inaccuracies  in  the  calculated  scattering  behavior  are  produced  in  the 

I 

(approximate)  evaluation  of  the  physical  optics  integrals.  In  the  p.-esent  case, 
however,  this  approximate  evaluation  yields  resmu  which  are  quite  acceptable 
for  many  purposes.  The  qualitative  agreement  between  the  compmient  echoing 
areas  is  good  (see  figure  5-3)  and  indeed,  the  approximate  values  are  in  error  by 
no  more  than  10  per  cent  for  a  ^here  of  radius  5)(/v  providing  the  receiver  lies  in 
the  half-8par«  cmitainiag  the  incident  field.  Even  if  the  bistatic  angle  exceeds  90° 
the  eiT' -c  in  using  the  optics  formulne  remain  small,  and  for  a  =  a  permitted 
error  of  16  per  cent  would  extend  their  validity  to  cases  of  scattering  throu^  an^es 
as  large  as  120°  . 

is  view  of  this  agreement  there  seems  every  reason  for  ptitting  forward  tiie 
physical  optics  scattering  cross  section  for  use  in  practical  calculatimis  invdlvii.g 
q[>heres  of  radius  greater  than  3\/2,  Providing  the  receiver  is  directed  at  the 
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iiAuminated  portion  of  the  sphere,  the  component  scattering  cross  sectio;^  for  a 

linearly  polarized  incident  plane  wave  are 

2  2. 

Oq  =  sra  cos  p 

2  ?. 

CTjj  =  Ta  sin '  0 

where  6  and  0  are  polar  coordinates  defined  with  reference  to  the  directitms  of 
incidence  and  of  the  Incident  magnetic  vector.  These  results  are  sufficient  to 
define  the  apparent  cross  section  applicable  to  any  receiving  system  and  for  any 
type  of  incident  polarization. 
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